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Abstract Starting from yRON’s vapour pressure equation, separate generalized correla- 


tions are developed for the gas and liquid phose contributions appearing in this equation. The 


JOL. fine) equation obtained by combining these correlations is : 
12 L $58 log p, OO244 (p, 1) 
960 T 


T Ty Tr 
| (1/277) dT, (w 
i 
\ rigorous graphical solution of this equation is presented which enables estimations of the heat 
of vaporization of a substance from a knowledge of its critical tempcrature, critical pressure 
and vepour pressure dvta. The correlation hes been tested for a vericty of substances up to 
the critical point and found to be very accurate, porticularly for orgenic compounds and elemen- 
tary (or simple) molecules for which the averege error in the use of this correlation is less than 
2 per cent. 


Résumé—Partant de léquation de CLarryron, les auteurs développent et généralisent des 
relations individuelles relatives a la contribution des phases gaz et liquide que cette équation 
met en jeu. L’équation finale obtenue a partir de ces relations est : 
5 
4.58 log p, 0.0244 1) 
Ty Ty 
i t 
Une solution graphique rigoureuse de cette équation permet destimer la chaleur de yvaporisa- 
tion dune substance a partir de sa température critique, de sa pression critique et de sa pression 
de vapeur, 
La relation a été vérifiée pour plusieurs substances au dessus du point critique, et clle s'est 
révélée trés précise, en particulier pour les corps organiques et les molecules simples, pour lesquels 
lerreur moyenne est inférieure a 2 pour cent. 


Zusammenfassung —\usgehend von der Dampfdruckgleichung nach Clapeyron werden 
yetrennte verallgemeinerte Beziehungen fiir die Anteile der Ges- und Flissigkeitsphase in dieser 
Glei hung entwickelt. Diese Bezichungen werden zu folgender Schlussgleichung zusemmen- 
pesetzt : 

L $58 log p, (p, 1) 

T, 

(1/2T,*) dT, (w/T,) dT, 
i i 
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Kine genoue graphische Lasung dieser Gleichung wird mitgeteilt, welche die Abschitzung der 
Verdampfungswarme einer Substanz aus der Kenntnis ihrer kritischen Temperatur, ihres 
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kritixchen Druckes und der Dempfdruckdaten gestattet. Die Bezichung + urde an zahlreichen 


Substanzen bis zum kritischen Punkt gepriift und fiir sehr genau befunden, insbesondere ftir 


organische Verbindungen und clementare (oder 


Fehler kleiner als 2 Prozent ist 


SuRVEY AND OnsectTivi 


M ANY empirical methods have been propose ad im 
the past for the estimation of latent heats of 
vaporization at different temperatures and pres 


sures. The familiar Trowron’s rule and its 


earlier modifications [1, 2] are of limited use im 
process ca ulations: the modifica 
tion [3] of the Trowuron rule has merited some 
attestion [4, 5] but has not vic Ided any generalized 
relatio. ‘hip. A simple equation was proposed by 
KisTiaki wsky [6], but it is applicable to non 
polar hq? ids at their normal boiling points only. 
Another noteworthy relationship is the Antone 
equatic a of and this equa- 
tion ‘avolves two constants and. im addition 
requin s vapour pressure and critical point data 
for its . se. 

The m. st useful, if not the most accurate, 
method ava ‘able is perhaps that of WATSON 
[8, 9}. Based o. the observation that, when the 
TrovTon ratio is , lotted agaist red wed tem 
perature, curves ot tt sar shape ire obtained 
for all substances, whether polar or non-polar, 


Warsow obtained the general « mipiri il equation, 


which can be used to obtain the heat of vaporiza 


tion at any given temperature and pressure 


provided the value at any one temperature and 
pressure (usually the normal boiling pomt) ts 
known together w’ ‘ that of « ritical temperature. 


This method has the advantage that vapour 


pressure data are not required, but large errors 
may be encountered in the vicinity of the critical 
point. 

The well known method of reference substance 
plots proposed by Orumen [10] was later improved 
by Gorpon [11) who replaced the equal-tempera 
ture basis of these plots by the more general 
equal-reduced temperature basis. Keceently 
Cur Li and Cansan [12], utilizing the observation 
of LypDERSEN et al. [13] that the difference between 


cinfache) Molekille, fiir welche der mittlere 


gas and liquid phase compressibilities ts a unique 
function of reduced pressure only, proposed a 
reference substance method based on equal 
reduced pressures and reported good agreement 
between calculated and observed heats of 
vaporization. 

\ ralized procedure, employing reduced 
co-ordinates, was developed by Meissner [2, 14). 
On the assumption that the relationship, 

log p 
is rigidly valid up to the critical point, Meissner 


derived the equation, 
low P, 


As pointed out by Warsow [9], the use of equa- 
tion (3) becomes somewhat unsound at low 
values of T,. 

The object of the pre sent paper ts to offer a 
ric orous, alternative »rocedure using generalized 
co-ordinates for an accurate estimation of heats 
of vaporization up to he critical point, and to 
present graphical soutions of the equations 


developed, 


TurorericaAL DrveLorMEN' 
The Clapey ron equation in its original form ts 
written as: 
dp 


dp 


dT dd 

The latent heat may therefore be visualized as 
consisting of two contributions, one from the 
gas phase and the other from the liquid phase. 
Denoting these contributions by L, and Ly, 


respectively, we obtain 


L=L, (6) 
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L, and L, may now be separately evaluated and 
combined to obtain L. 


Evaluation of L, 
From the following two equations, 
dp 
dT 
and RT (8) 


(7) 


L, =vT 


equation (9) below may be obtained by substitut- 
ing reduced temperature and pressure (7', and 


p,) for T and p and separating the variables : 


L, 
- — dT dp : 
RT, (9) 


r 


Equation (9) can be integrated if the limits of 
integration are fixed. This is easily done if it is 
remembered that the latent heat of vaporization 
of a liquid is zero at the critical point ; the limits 
may therefore be written 
and T, 


Thus, assuming as a first 
approximation, an average constant value of 
L, over a given range of T,, values, 
2.303 R log P, 
Tr 
2 "yr 
1 27 dT, 


(10) 


Ty 
When the term (dT, =T7,*) is integrated on the 


i 
further assumption that zc is constant at an average 


value over a particular range of temperatures and 


pressures, we obtain 


L, 2.303 log p, 
(1 1/T,) 


(11) 


which is similar to the Meissner equation for L. 

But it is evident from the assumptions made that 

equation (11) is unsound, Therefore, in order to 

arrive at fairly accurate values of L,, it is necessary 
Ty 

to evaluate the term ( dT, =T2 by graphical 
i 

integration. This procedure is elaborated in a 


later section. 


Evaluation of 


The liquid phase contribution is given by 
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dp 


(12) 


If v, can be expressed in terms of 7, and p,, 
equation (12) can be generalized. For this 
purpose, Warson’s postulation of the liquid 
expansion factor [9], defined by equation (13), 
may be used : 

p. M v, (15) 
Watson constructed a chart for w as a function 
of T, and p, using the experimental values of 
iso-pentane, normal pentane and propane. 
Bresacn and Tuopos [15] have recently con- 
structed another chart based on the experimental 
values of nitrogen. The two charts show devia- 
tions of the order of 10-15 per cent. This devia- 
tion is explained by the fact that the values of 


the critical ratio, 


x, = (14) 
RT. 

are different for the liquids used in the construc- 
tion of these charts; the average value of the 
three liquids employed by Warson is 0-271 and 
that of nitrogen is 0-291. It has been suggested 
by Bresacu and THopos [15], Lypersen ef al, 
[13] and others [16, 17] that in addition to T, 
and p,, z, should also be employed as a correlating 
parameter. Consequently, neither of the two 
available w-charts can be considered as truly 
general. This may now be examined with respect 
to the present study. The errors introduced by 
neglecting variations in may be obviated by 
the observation that the product (w v,) is constant 


for a given liquid : 
wi 1 k { 15) 


This may be combined with the further observa- 
tion that if v, at a reference state is known then 
v, at any other state may be calculated with 
accuracy by the equation, 

(16) 


w 


using a chart of w as a function of T, and p, 


prepared for any liquid. This is done by first 
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calculating ( Up )reg by use of the w-chart and 
then obtaining v, by reading the value of w at 
any desired temperature and pressure from the 
chart. In this way the errors resulting in the 
direct use of any available w-chart for liquids 
having different values of 2, are reduced to less 


than 5 per cent. 


As mentioned above. the i. of equation (16) 


is dependent on the value at a reference state 
being known. The choice of the critical point as 
the reference state generalizes this 


At the 


nearly 


equation, 
critical point, the value of w is very 
constant at 0-044: substitution of this 


value in equation (15) gives 
k (17) 


Combination of equations (14), (15) and (17) and 


rearrangement results in equation (18 


r 
v, = 0-044 z R | (18) 

aw 
The actual fixation of the value of = to be 


employed in equation (18) can be done by taking 
the best average for a variety of liquids, This 
value has been found [18] to be 0-276 which is 
very close to the average value for the liquids 
employed by Warsox. Substituting 0.276 
in equation (18), combining this with equation 
(12) and placing the limits of integration as 
before, one obtains by use of yveneralized co- 


ordinates : 


1 
aT, 00122 R(p,—1 (19) 
The term |(L, @ T,)d 7, may now be integrated 
; 
using the w-chart by assuming L, to be constant 


For values of 7. 0.6 


and for low pressures, w has been expressed as a 


as a first approximation. 


straight line function of 7’, ; 


O-1745 0-0838 7) (20) 


Substituting this in equation (19) and integrating, 


L, — 1) 


0.0838 (7° 1) 


21 
0-408 log T, 
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For higher values of 7, equation (20) is no 


Tr 
longer applicable and the term (L, T,) dT, 
J 


i 
should be graphically integrated as discussed in 


the following section. 


GRAPHIC SOLUTION 


The generalized compressibility factor chart 
constructed by Hovucen and Warson [19] may 


be used for the integration of 


} (1 27) dT, of equavion (10), 
i 

assumed 7°. values up to unity, graphical integra- 


graphical 
For a series of 
tions were performed for p, values from 0-03 to 


O-975. 0.03, the values of 2z 
between 0-98 and 


being 


For P, 
approximately 10, 2 Was 
assigned an average constant value for a viven 
p,, and L, was calculated directly. The chart 
for L, T,. 


presented in Fig. 1. 


showing only the main values of p,, is 
It may be pointed out that 
graphicel integrations performed as above would 
actually 


“ive a mean value L, and not L,, L, 


being defined by 


oT? / 


. r . 


\ few 


ranves at 


integrations over different temperature 


chosen constant values of p, have 
shown that L, is practically independent of as 
so that L, ~L,. Thus Fig. 1 may be considered 
quite accurate. 

For evaluating L, T,, the analytical solution of 
equation (19) presented in equation (21) was 
combined with graphical integration of the term 

\(L, a T,) dT, using Warson’s chart in order to 


obtain values of L, 7, up to the critical point ; 


these are presented in Fig. 2. It has been 
observed, as in the case of L,, that the mean 
values, L,, are equal to L,. Fig. 2 shows the 
values of L, T, over the entire range of interest : 
the nature of the curves suggests that two 
different graphs are required for accurate reading 
of the values, but in the figure these are combined 


by a break of scale and only the main parameters 
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hic. 1. Generalized correlation for L, ao 


are drawn. For values of 7, > 0-7 the constribu- while for values less than 0-7 this contribution is 


tion due to the liquid phase is seen to be sizeable, very small. 


| 
30 
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4 Pia. 3. Generalized correlation for heat of vaporization T KO OL cal g mole 
The final chart for estimating L. 7, is presented 
in Fie. 3. This was constructed by enlarging 
various sections of Figs. 2 and 3 and combining -/ 
the values read from them according to equation * 
(6). Curves have been drawn for a large number ho 
of p, values to enable accurate interpolation. 
While the limits of the parameters in Fig. 1 | 
have been arbitrarily tixed, those in Figs. 2 and 3 4 a’ 
A 
: have been fixed on the basis of the expansion = 4000 | rs | 
factor chart and available liquid compressibility 
data. 2000 
Accuracy of correlation 
0 4000 
The correlation developed in this paper has o0C 8000 
L(experimenta!) 
heen subjected to an exhaustive test. Values 
read from Fig. 3 were compared with ¢ xy rimental Fic. 4. Comparison of experimental and calculated 
results for a variety of substances which included heats of vaporization 
hvdrocarbons, organic acids, bases, esters, ethers, 0-009 0-990 T, 0-500. 0-990 
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halides, elementary molecules and a few inorganic NOTATION 


compounds, A total of about 100 comparisons 1, B — constants of vapour pressure equation 

(up to the critical point) were made and the k = constant of equation (15) 

correlation was found to be good to within about L. = heat of vaporization cals/g mole 
2 per cent on an average. Maximum deviations L’ = heat of vaporization at a reference condition of 
were encountered in the case of a few inorganic temperature 

substances which gave errors of the order of L = mean value of L defined by equation (22 
10-15 per cent. In these cases the other available VM = molecular weight 

correlations gave errors of a higher magnitude. p = vapour pressure 

Fig. 4 shows a plot of calculated and experimental R = gas constant 

results ; to avoid heavy clustering of points, only T — absolute temperature K 
about sixty sets of data are plotted covering the T’ — reference temperature 


entire range up to the critical point, compressibility factor 


On the basis of a survey of all available methods : angles welutitn 


for estimating heats of vaporization, it is concluded 


w liquid expansion factor 


that where vapour pressure data are not available 


Watson's correlation [9] appears to be the best Subscripts 
apy 


known empirical procedure, but where vapour l gas phase contribution to heat of vaporization 
pressure data are available the rigorous develop- 2 liquid phase contribution to heat of vaporization 
ment presented in this paper should give very ‘ critical state 

accurate results, particularly for organic com- l liquid state 

pounds and elementary (or simple) molecules. r = reduced state \ 
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An investigation of the viscous and inertial coefficients for the flow 
of gases through porous sintered metals with high pressure gradients 


D. B. Greenperc* and Wercer 


Department of Chemical Engineering, The Johns Hopkins University, Baltimore 18, Marvland 


( Received 12 May 1850) 


Abstract — This paper deals with the results of an investigation of the steady-state flow of gas 
through porous sintered metals under high pressure gradients. The purpose of the study was to 
determine the behaviour of the viscous and inertial flow coeflicients in such media as a function 
of pressure and temperacure Pore size distributions of the various cores tested were obtained 
by means of mercury intrusion measurements. Mean cross-section diameters which were derived 
from this data correlated well with experimentally obtained permeabilities. The porosities of 


all six of the cores which were used were within the 30-40 per cent porosity range. 


It is found that the permeability (the reciprocal of the viscous flow coeflicient), while not a 
function of pressure, decreases significantly as the temperature of a sintered porous medium is 
increased. This is attributed to small changes in the surface structure of the particles which 
constitute the medium. The inertial coeflicient for any given core is found to be a constant within 


the pressure and temperature ranges which were covered by the present experiments. 


Résumé —-Ce mémoire traite des résultats dune étude sur lécoulement permanent d'un gaz a 
travers les pores de metauyx frittes sous de hauts gradients de pression. L’objet de cette étude 
etait de determiner le comportement des coeflicients de viscosité et dinertie de Técoulement 
dans ce milieu en fonction de la pression et de la température. Les distributions de la dimension 
des pores des differents corps examines étaient obtenus au moyen de mesures par injection de 
merceure Les diameétres movens des sections transversales qui sont obtenus de ces données se 
relient bien avec les permeéabilités obtenues expérimentalement. Les porosités des six corps 


utilisés se situent dans Vintervalle de 30-40 pour cent de porosite. 


Zusammenfassung Die Arbeit beschiftigt sich mit Untersuchungsergebnissen der gleichfér- 
migen Stromung von Gas durch pordse Sintermetalle unter hohen Druckgradienten. Zweeck der 
Untersuchung war, das Verhalten von Viskositits- und Trigheitskoeflizienten in derartiven 
Stoffn ats Funktion von Druck und Temperatur zu ermitteln. Die Porengréssenvereteilungen der 
verschiedenen gepriften Kerne wurden durch Infiltrationsmessungen von Quecksilber ermittelt 
Die so erhaltenen, mittleren Querschnittsdurchmesser stimmten gut mit den experimentell 
ermittelten Permeabilititen tiberein. Die Porosititen aller sechs verwendeten Priifkérper lagen 


innerhall des 30 bis 40 prosent Porosititsbereiches. 


Es wird festgestellt, dass die Permeabilitét (reziproker Wert des Viskosititskoeflizienten), 
die nicht vom Druck abhiingt, betriichtlich mit steigender Temperatur des pordsen Sinterwer- 
kstoffes abnimmt. Das wird den kleinen Anderungen der Oberflichenstruktur der Teilchen des 
SinterkOrpers zugeschrieben. [es wird gefunden, dass der Triigheitskoeflizient fiir jeden Prifkorper 
innerhalb des Temperatur- und Druckbereiches, tiber den sichd ie Untersuchung erstreckte, 


konstant ist. 


*Present Address: U.S. Naval Academy, Annapolis, Maryland. 
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INTRODUCTION 


Tue Neep for study of the flow of fluids through 
porous media has become apparent in many fields 
of engineering and science in recent years. An 
understanding of this phenomenon is important 
in the study of ground water hydrology, soil 
mechanics and natural gas reservoirs of the 
petroleum industry. Much of the early work has 
been compiled by Muska, first in 1937 [1], and 
again in 1949 [2], and by Hurserr [3], in 1940. 
Related work in absorption, extraction, distilla- 
tion and filtration has also yielded important 
results in this field. Current interest in the flow 
of fluids through porous materials is centered 
about applications in aerodynamic boundary 
laver control and transpiration cooling in the area 
of ballistics and rocketry |. 5). 

The steady-state flow of fluids through porous 
media has been studied extensively from both an 
experimental and = analytical standpoint for 
relatively low pressure gradients. GRUNBERG and 
Nissan [6] investigated viscous effects in gases at 
low tlow rates. The deviation of flow in consoli- 
dated media from the characterized by a purely 
viscous motion has been studied by Green and 
Duwrz [7]. Correlations for both viscous and 
inertial effects in unconsolidated media were 
presented by Browne et al. [8]. 

The conclusion to be drawn from these papers is 
that most experimental results for one dimensional 
flow through porous materials (both consolidated 
and unconsolidated) may be correlated by an 


equation of the type 


+ Bp V™ (1) 
da 
where « and 8 are so-called * constants,” and n 
is some number between 1 and 2. The parameter 
x is the inverse of the famous Darcy constant, 
k, [9]. The first term on the right-hand side of 
equation (1) represents the portion of the pressure 
eradient due to viscous effects, while the second 
term represents the contribution by inertial 
effects. The exact nature of the inertial con- 


tribution is not well understood, It is unlikely 
that true turbulence exists in a_ fine-grained 
porous material, The majority of the inertial 
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effects are probably due to expansion and con- 
traction of the fluid, rapidly changing direction 
of the flow path, ete. 

The “ constant” 8, as well as x, is, of course, 
dependent on the geometry of the medium. While 
it is impossible to describe this geometry in its 
microscopic detail, except for idealized media, 
quantities such as porosity, specific surface and 
pore size distributions — as well as the structure 
of the component particles—can give a gross 
indication of the structural variation among 
various porous materials. The structure may 
change to a certain extent as the ambient con- 
ditions are varied. Thus, changes in pressure or 
temperature might be expected to affect the 
values of « and,or 8. Little mention has been 
made in the literature of these effects. 

Cornet and Karz [10] have presented correla- 
tions of data on the flow of gases through con- 
solidated porous material using equation (1). 
Their results indicate that n = 2 for these types 
of materials. Under steady-state conditions the 
mass velocity of gas (for one dimensional flow) 
through a porous material must be constant. 
If we let pV = G, equation (1) becomes 


dp 


—p— + pG (2) 


da 
Using the perfect gas law with a compressibility 
coetlicient, we obtain 


M dp d(p*) 
P da da (3) 
Integrating and rearranging equation (3), we 
obtain, 
A (p*) Mg G 
2L (zRT) 


which is the relationship CorNeELL and Karz 
used in analysing their data. If the quantity on 
the left-hand side of the equation is plotted 
versus G/y for a given run, the slope will yield 8, 
while « will be obtained from the intercept. 
The Reynolds number of a flow system is 
defined as the ratio of inertial to viscous forces. 


Since the ratio of 8 x represents the ratio of 


the inertial to the viscous coefficient, a modified 
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Reynolds number for porous media may be flow equations for porous media leaves open the 


detined as question whether the “constants” are really 


(")(") . constant and vary only as the structure of the 


fe medium varies: or whether they themselves are 
By inspection, equation (4) may be rearranged dependent to an extent on the pressure and flow 


to yield : rate obtaining in the medium. Recent interest 


in the high pressure region, both for steady-stat« 

A (p*) Mg ‘ and transient flow in porous materials, has 

(2RT) a) (Gp demonstrated a need for investigations of the 
applicability of the relationships for flow through 


The left-hand member of equation (6) ts a ratio 


ofl dissipative to the rtial forces and the retore, porous media white h have been le ve loped mainly 


represents a friction factor. Thus, equation (6) 0" the basis of low pressure data. The constancy 


of the “ constants ” especially must be checked. 


becomes 
It was the purpose of the present study to 


obtain ¢« xperin ntal evidence on this point overa 


yressure range up to ZOOO p.s.i. 
and has been suggested by Green and Diwerz I 


17 to repre sentative tmnetion factor corre la Expt! RIMENTAIL 


tior media, Tl riction factor 
on for porous mec Phe frictio cor Six porous sintered metal cores, covering a wide range of 


repre sented by 2 Vin. in the Darcy (or laminar metal particle size and permeability, were selected for this 


re“ion and by 2 in the mertial region investigation. These samples were chosen after a series of 


Due to lack of knowledge of the exact fin preliminary tests had determined the characteristics of 
structure of most porous media und the metallic media which were necessary for obtaining re- 
producible ind reliable data The powder size and sinter 
inadequacy ol present analytical methods for 

ing conditions for some of the cores were selected in order 
describing fluid flow through such complex to obtain the high pressure gradients which were of primary 
structures it is necessary to re-determine the interest in the present work. The cores which were used 
validity of EXPTESSIONS St h as equation t) when for the reported runs and their physical characteristics 


any of the parameters are changed drastically are listed in Table | 
; The particles (except for the brass core) were approxi 


from those which previous xperimental > 
mately spherical in shape. Pore size distribution measure 
data was obtained. ments were performed on samples eut from the cores by 


The unrigorous nature of the derivation of the means of a mercury intrusion porosimeter [15] Fig. 1 


Table : Physical properties of porous cores 


Porosity (22 
pore 
diameter Krom 
Particle siz (a) density pore sim 
Supplies ” ore range (from pore asurements mreasurcments 


measurcments) 


Frankfort - Bronze 
Arsenal Brass 


Kulite Nickel 


Tungsten Co Copper 


Svivania Nickel 


Electric Corp T Tungsten 


10 


12 

|_| 


Total volume | 

22°C 0.052 at 22°C 

0-04 e 42°C 0-050 ot 42°C | 

5 | 
0-03 | 


° 
O 


Pore volume, 


Pore diameter, 


1. Pore size distributions of the S core, 
Total volume 


o 22% Oos2 at Zz ¢ 


shows distributions of the 8S. N-—1 core at two temperatures. 
Similar results were obtained for the other cores. The 
reason for measurements at two temperatures will become 
apparent later. The data were a tually reported in terms 
of fraction of pore volume over a given size range. The 
points in Fig. 1 represent the midpoints of these ranges 
Mean pore diameters were determined from these curves. 
The method, based upon the principle of capillary forces, 
vields an approximate picture of the internal structure of 
a porous material. However, it is important to understand 
the limitations of the method. In a material such as 
sintered metal there are many spaces or “ cells which are 
connected to the rest of the porous system only by passages 
which are of smaller average cross-section than the cells 
themselves. A method based upon the increase in pressure 
required to force a liquid, such as mercury, into smaller 
and smaller pores will include the volume of such cells 
with the pore size range which contains the diameter of the 
largest passage connecting the cell and the rest of the flow 
system. The results are, therefore, essentially measure- 
ments of the “ constriction " diameter distributions rather 
than true average pore size diameter The importance of 
these contributions will be touched upon further in the 
discussion 

The other limitation — one which is not important in 
the present work —is the fact that one can only measure 
the pore space which is connected by tlow passage to the 
rest of the system. Completely enclosed cells will not be 
included. It can be seen in Table 1 that the void fractions 
obtained from the pore size distribution measurements are 
usually below the corresponding values determined from 
density measurements. It was determined that no 
appreciable structure of the cores used in the recent 


experiments lay outside the 2000-1 limits of the 
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measurements. The difference in densities must, therefore, 
be due to enclosed cell structure. An anomalous result 
was obtained with the S-T-1 core. At present, there is 
no satisfactory explanation for this case. 

A potting method was developed for encasing the 
porous cores in steel sleeves using an epoxy resin formula- 
tion similar to those used by manufacturers of electrical 
components. Because of the extreme wettability of this 
resin, the cores (all except the S-T-1 core which was 
‘sweated into its sleeve) were sealed in an aluminium 
putty-like resin, to prevent the epoxy resin from seeping 
into the pore spaces of the sample. 

After the samples were mounted in the sleeves, both core 
faces were given finish cuts on the lathe to make the ends 
flat and parallel. This operation tended to smear the 
porous metallic surface necessitating a me tallographic 
treatment to restore the end grains as closely as possible 
to being representative of any cross-section within the 
sample. 

This metallographic process included a wet sanding using 
300 W and 600 W silicon carbide polishing papers followed 
by a final treatment on a series of alumina polishing wheels 
to bring out the grain boundaries, preparatory to opening 
the pore surfaces. 

An electrylotic polishing method was used to remove 
the remaining smeared metal from the surface. It was 
found that the ele« trolyti« technique was far superior to 
the more common chemical etching methods for this work. 
Fig. 2 is a photomicrograph ( * 250) of the surface of core 
F-B2-12 (Bronze) after such treatment, while Fig. 3 
shows the surface of Core F- Br4 32 after similar treatment. 
Details of the technique are reported in references [11, 12 

The experiune ntal apparatus consisted of an hydraulic 
svstem designed to maintain a constant source of high 
pressure nitrogen gas for the flow studies, a gas reservoir 
and a suitable means for controlling and measuring 
pressure, temperature and gas flow rate. Fig. 4 is a flow 
diagram of the equipment 

The pump was a dual feed positive displacement pump, 
with a continuously varying stroke control to regulate the 
hvdraulic fluid output up to a maximum of 4-9 gal hr for 
each feed, at 2500 p.s.i. The hydraulic system functioned 
continuously during a run supplying oil to the gas reservoir 
at a rate slightly above that of gas withdrawal. The 
excess oil was bypassed back to the storage tank. This 
arrangement insured positive control while damping 
pressure surges inherent in positive displacement pumps. 
\ solenoid-pneumatic piston type control valve metered 
excess fluid back to the storage tank. The solenoid was 
controlled by a pressure-sensitive electric contact gauge. 

Adequate control over the gas flow system was obtained 
by a gas-loaded diaphgram valve. Fluctuations in gus 
flow were observed to be generally less than about 2 per 
cent. The sleeve containing the porous core was fixed 
hetween two end blocks which contained the connections 
for gas inlet and outlet and the taps for pressure measure- 


ments at both ends of the sample. A satisfactory pressure 
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pressure gauges gas lines 


oil lines 


® valves 


flow meter 


nanometer 


systematic investigation of the temperature dependence 
of permeability 

maintained 
\ 
pair of 200 w flexible heating tapes, one wrapped around the 
the 


other around the c« 


Temperature control and regulation was 


by means of a water bath and flexible heating tape. 


line from bath to the porous core assembly and the 
itself, were employed to facilitate 
the test 
was adjusted individually until the inlet and 


two iron constantan thermocouples 


re 


temperature control in device. The current in 


each 
outlet 


temperatures agreed within Temperatures 


were measured 
placed about 3 in. upstream and downstream from the ends 
of the sample 

dual 
this 


meter was better than 0-5 per cent of the full scale reading, 


Gas tlow rate measurements were made with a 


clement orifice meter Accuracy obtained with 


the orilice was designed to give a linear flow rate 


sites 


pressure drop relationship. Calibrations were made 
infrequent intervals with a wet test meter, 
All 


safety barricade of steel boiler plate 


high pressure equipment was housed behind a 


All valves, measuring 


instruments and controls were mounted on the outside 


so as to facilitate operation from the front of the shield. 


The high pressure gas reservoir was fed from a four- 
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cylinder gas manifold mounted outside of the barricade. 


Some runs were attempted at lower temperatures by 


packing crushed ice around the sample and filling the water 


bath with an ice-brine mixture. At high flow rates. 


however, it was rather diflicult to maintain a constant 


temperature by this method ; consequently, data from 


only one sample has been reported at a temperature 


below ambient conditions. 
Further details of the experimental apparatus and 


method are available in ref. {12). 
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A series of four to six runs at different tempera- - — . 
tures was made for each of the six porous cores, > 256° 
The temperature range was between 15 °C and 
60°C. Pressure gradients ranged from 790 p.s.i. in. 
for the S—-T-1 core to about 25 p.s.L./in. for the 
K-C--1 core ; these corresponded to a volumetric hic. 6. Cornene and Karz plot for K C1 core 
flow rate of approximately 100 ft® hr, which 
represented the maximum obtainable with the os dining 
flow measuring apparatus. Fig. 5 is presented 
as an example of the type of data obtained. In the majority of cases, the data were easily 
This is the data for the K-C—-1 core. The complete fitted by eye to a straight line, as predicted by 
data is available in ref. [11, 12}. equation (4). However, for certain cores, there 


was an upward deviation of the points in the low 


ra flow rate region at elevated temperatures. For 
fh. these cases, the higher flow data points were used 


a to define the curves. The reasons for this are 


| Ia explained in the Discussion. 

KEE The values of x (and the corresponding permea- 
go bility, &) and 8 obtained for the six cores are 


wresented in Table 2 as functions of temperature. 
| 


Calculated permeabilities range from about 


+ millidarcies (mdar) for core S—T—-1, having the 


sc.fh /in 


Fic. 5. Pressure drop vs. flow rate curve for KC. 


Temperature, 


core, 30} 
@ 57-6 @ 447 
20} 


| 
| 


Values of x and 8 were determined according 
to equation (4). As an example of the Corners 10 2 4610 2 4610 2 460 2 4610 
and Karz plots which were made, the plot Permeability,  mdor 
corresponding to the data of Fig. 5 is given in Fic. 7. Vermeability vs. temperature for porous 


Fig 6 cores 
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Table 2. Values of x. 8B andk 


147710) } (mdar)* 


lise 

835 


*1 1-06 ft? 


smallest mean pore dimension, to 3000 mdar for Ak k 
core K-C-1, which had the largest mean pore AT 
size value. Fig. 7 is a graph of permeability 

plotted against temperature for all of the samples where Ak k k (60 ©) k (20 C) k (20 °C) 


studied. These curves were fitted to the data by 


per cent change ink per © 


the usual least squares technique, assuming that and AT i 

the relationship was linear over the range of 

temperatures investigated. In all cases except Fig. 8 is a representation of the data according 
core 5 T-1, permeability decreased with increas- to equation (7). The values of x and 8 presented 
ing temperature. The percentage change of in Table 2 were used for the calculation of the 
permeability with temperature has also been Reynolds numbers. The line is the (, calculated 
tabulated from the graph in Table 3 as : from equation (7). 
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Friction factors vs. Reynolds number 
32 
S-N-1 


Table 3. Thermal coefficient of permeability Because of the Joute—Tuomson cooling effect, 
the gases were heated as they expanded in the 


Porosity (From 


Core density measurements ) \T cores so that the operation would remain essen- 


(%) ec-ly tially isothermal. This required that the tempera- 
ture of the solid be somewhat higher than that of 
the gas. In the flow regions where the heat 
transfer was reasonable eflicient, Le., where 
appreciable inertial effects existed, the tempera- 
ture gradient could be low and the required heat 
would still be transferred. At low tlow rates, 


however, the major amount of heat is transferred 


to the bulk of the gas by conduction and this 
Discussion means that the temperature gradient must 
It was mentioned previously that in some cases increase to maintain an adequate heat supply 
the experimental data deviated from straight to the gas. The temperature of the porous 
line behaviour when plotted according to the medium in this case would then have to be 
method of Cornett. and Karz. This occurred considerably above the temperature of the gas 
especially for the S-N-1, K-N1-2 and K-C-1 (which was the temperature which was actually 
cores in the low flow rate region at elevated measured at the inlet and outlet to the core 
temperatures. The effect can be attributed to the — since the actual core temperature would have 
heat transfer which took place between the metal — been diflicult to obtain). This conclusion ts borne 
and the gas in the cores. out by the fact that temperature control in this 


Vx | 
| 
VOL. | | 
| 
s 
K 
S-T-1 
S-N-1 : 
K Ni-2 
KC 
F-B2-12 | 
F_-Br4—32 ay 
15 


D. B. 


region was considerably more diflicult than in 
other regions, 

The curves which would be obtained by drawing 
a line through the flow rate points in the higher 
K-N1-2 and 


K-C-1 cores would, therefore, not be expected 


temperature region for the S-N-1, 
to represent isothermal core conditions. There- 
fore, in these cores, in view of the lack of data on 
actual core temperature, it was considered more 
reliable to take the data at the higher flow rates 
as representative of the conditions of the cores at 
temperatures close to that of the gas. These data 
vave a straight line behaviour which was consistent 


with the data at lower temperatures. This linear 


behaviour indicates k and 8 are independent of 


flow rate or pressure for the materials under 
consideration, 

Examination of Table 2 reveals the behaviour 
of k and 8 
in all except one case, independent of temperature 


The B2 


there is an abrupt 


as a function of temperature. 8 ts, 


within experimental error. 12 data is 
anomalous in that above 27 
increase in 8. This may have been due to a 
structural change with pressure, or the penetra- 
tion of some unknown obstruction into the core, 
Despite these apparent discrepancies, however, the 
the that the 


inertial coetlicient 8 does not vary appreciably 


weight of evidence seems to be 
with temperature for any given core. 
the the 


temperature dependence of the permeability Is 


In view of invariance of 8, strong 


remarkable. This is an effect which had not been 
previously reported in the literature and may 
explain some of the “ scatter” in permeability 
data for sintered metals which has been reported. 

In all except one of the reported cases, & 
decreases with temperature. To arrive at a 
reasonable explanation for this effect, we must 
consider the physical significance of & and 8. 
In some complicated fashion the permeability 


is a measure of the free cross-sectional area and 


length of flow path in the porous media. (These 
quantities, of course, consist of integrated averages 
for a macroscopic section of material). In other 


words, the reciprocal of &, x, is a viscous resistance 


coeflicient. Let us restrict ourselves for the 
moment to one hypothetical “path” the 
porous material, The fluid in the “path” (a 
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the 


passage is not enclosed) will encounter successively 


quantity which varies, of course, becaus 


portions which have relatively large  cross- 
sections followed by narrower passages near the 
points of contact of the spherical particles which 
constitute the porous sintered medium. Velocities 
will be higher at and near these “ constrictions,”’ 
and, therefore, the largest portions of the pressure 
drop along the path will occur there. The major 
the 


resistance coetlicient will, therefore, also orivinate 


contribution to overall measured viscous 
at these constrictions, 

As was mentioned previously, the pore size 
distribution measurements are primarily measure- 
We can calculate 


an average pore size from these data by taking the 


ments of these constrictions. 


it, and 
The 


square root of the sum, weighted by the appro- 


midpoint of each size range, squaring 


multiplying it by the appropriate volume, 


priate total volume. represents an average pore 
diameter. It is a rather rough average, because 
of the approximations involved in taking the 
midpoint of each range as the representative 
Nevertheless, 


the magnitude of the 


diameter. should ive some 


indication of * constric 


tions " we are dealing with. It would be expected 


that the permeability coeflicient ts related to this 


average diameter — at least for simular types of 


media, That this is, indeed, so is evident from 


Table 4. 


Table 4. 


Comparison of permealilities and mean 


pore diameter at approximately 22 
Cor \ (a) 
S-T-1 
K-Ni-2 
F Bz lz 3 
K 


While it is not possible to derive a simple 
functional relationship between & and y D*, it 
can be seen that the first five cores in the above 
table follow the same progression im kas im \ 1%. 


The F- Bré 32 core, being composed of irregularly 
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shaped particles rather than spheres, would 
rot be expected to fall in the same progression, 

The factor 8 is a quantity related to the relative 
amount of resistance to tlow existing in the porous 
medium due to inertial effects. In such a medium, 
these effects can arise due to several phenomena. 
The continuous expansion and contraction of the 
Inertial 


resistance is also caused by the many curvatures 


vas streams would ce rtaimly be one, 


in flow path as the gas wends its way through 
the material. Another important source of these 
effects ts surface roughness, just as an increase in 
the roughness of a pipe wall causes imereased 
pressure drop for tlow through the pipe. The 
amount of roughness ina sintered porous medium 
difficult: to define. It is probable 


ofl 


that a certam amount is) present on the 
particles before sintering. The sintering process 
causes small blisters and other surfac irregulari 
ties to artse in the surfaces, especially near the 
pomts of metal-to-metal contact. This is borne 
out by examimation of photomicrographs of the 
media. Since the velocities are largest near thes 
the amount of 


contact roughness in 


affect the total 


pots, 
these locations will markedly 
pressure drop through the medium. 

Since as discussed previously remains 
constant as the te mperature of a sintered Porous 
clement is varied, it ts apparent that its general 
structure remaims the same at least for a 


moderate temp rature range, That is. the 
contiguration of the tlow paths and the general 
character of the gas solid interface is not altered. 

Kxamination of the data reveals that the 
effect of the mertial and viscous terms on the 
pressure drop was of the same order of magnitude. 
The factor « is much larger than but, as can 
be seen in equation (1), it ts multiplied the 
viscosity of the gas a small quantity. 
Phe tive 


mvestigation which were mack 


porous cores used im the present 
from spherical 
powder particles all had porosities in the same 
range. It may. therefore, be assumed that. 
aside from surface irregularity effects, the internal 
structures of all tive samples were very similar, 
even after sintering. Examination of the pore 
size distribution curves bears this out. The 


regular relationship between increasing \ 1 and 
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increasing &, mentioned previously, is further 


proof, Excluding surface roughness effects, the 
3x ratio should, therefore, be the same for all of 
the cores under discussion, 

It is believed that the 


temperature ts mainly due to the presence of 


variation of k with 


these surface irregularities where the flow passages 
are narrow. Let us first. however. examine the 
feasibility of another possible ( xplanation of this 
effect. At first. it might be attributed to the 
normal thermal expansion of the metal particles 
which constitute the media. The thermal co- 
etheents of expansion forthe metals whose powders 
were used to fabricate the cores and sleeves used 
in the present study are tabulated in Table 5. 

Table 5.) Thermal coefficients of linear expansion 


for various metals 


Vetal Coefficient of linear rpansion we ¢ 
Copper 
Nickel 
2040 
Bronze 
Tungsten 
Stecl 


It is evident that the total thermal expansion of 
the metal particles over the temperature range 
used in the present study is negligible. Further- 
more, the coetlicients (except for tungsten) are 
all of the same order of magnitude. Therefore. 
it is apparent that the cause for the large tempera 
ture coellicent of k which varies much more 
widely between different samples than the thermal 
coctlicrent of « Xpansion, must lie somewhere else. 

As mentioned previously, it is believed that the 
root of the trouble can be traced to the surface 
irregularities existing at the points where the 
vreatest viscous resistance effects occur. Any 
slight change in the position or size of the irregu- 
larities can greatly affect the already relatively 
small opening available for the passage of the 
gas there. A small change in temperature of the 
medium, while causing no appreciable overall 
expansion of the particle, can cause a small shift 


plus small expansion effects in, 


in position of, 
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quant t with the 


tare 


look if the variation of thus 
materials 


lable 6 


some of 


thermal coetlicient of & for 


tested This data 


Is pore rite al iti 


While the 
due to the 


relationship i not ent consistent 


fact that two rent cs ol 


chill 

S N-1 and the F-B2-12 
cores, While their thermal permeability 
tendeney for \ 


chill 


are obtaimed for the 
oelherents 

to 
seen \s 


i! 


are the the 


mentioned previously ces 


cts. 


thre Cores composed ot spherical 


ribed mainly to differences surtace 


minute Ss in 


blisters 


onentation and or 


Apparently 


positions ol and other surtace out 


croppings due to temperature variations vreatly 
affect the 
This 


at positions in the cores wl 1 


permeability ot the int 
effect 


porous 
media would — be to be 
vreatest CNCTCIS 
nts at which 


reater the 


the largest influence on A 
the hiohest 


* the po 
velocities 


measured by thy 


will be thy 


ss effects. as para 


meter tha vreater temperature 


coetlicient of permeability. The overall change im 


porosity clin to this flect Is the his wis 
out by the distributor 


No 


temperature 


borne mcasure 


ments, significant changes in porosity. with 


wer noted for thre tor 


cores 
which measurements were made at 


he re 


» te ra 


tures, is evidence of 


amd 

size distribution with temperature (see Fig. 1), 
but at present it is not clear just how this change 
affects the phenomena which have been discussed. 
The 


determination of broad ranges of pore sizes could 


mereury mtrusion measurements, involy ng 


not be expected to pick up the fine variations 
which have been discussed, 

For the as before, the F-Brt-3z 
core would not be « xpected to fit into this scheme 
different 


used in its fabrication. 


reason 


because of the nature of the powder 
But the explanation for 
its large thermal permeability coctherent ts 
probably the same as set forth above. 

The 
of the S-T-1 
other coetlicrents, 
The 
into its steel sleeve rather than being “ potted,’ 


\s can bn 


Table 5 previous Vo pres nted, the linear coetlicient 


temperature coctlioent of rmeability 
to all the 


The reason for this is readily 


core opposite im sin 


apparent, tungsten core was “ sweated’ 


as were the other cores. seen from 


of thermal expansion of stecl ts approximately 
three 
20 the 


times that of tungsten. Assuming that at 


tungsten Was rfectly sealed the 


sleeve, any rise im temperature would cause 


greater expansion of the inside surface of the 


sleeve than of the outside surface of the tungsten. 
The 


was, therefore, 


increase of px rmeability of the S-T-1 core 
apparently caused by the creation 
of additional free cross-sectional area at = the 


tungsten-steel interface. 


CONCLUSIONS 


The present study has demonstrated that the 
inertial and viscous coeflicrents in the equation 
describing flow through porous media are indepen- 
dent of pressure at pressures up to 2000 p.s.i. for 
sintered metal « ompacts. It is, however, nec ssary 
to take into account variations of the permeability 
with te mperature even over narrow temperature 
ranges. A decrease in permeability as the tem- 
perature of the medium ts ratsed may also occur 
In Porous materials other than sintered metals. 
It is likely that the behaviour of porous surfaces 
such as are used in aerodynamic structures, and 
materials, such 


porous as solid propellants, are 


significantly affected by this phenomenon. 


It is also clear that more attention and study 


must be devoted to the internal surface structure 


Table 6. ¢ paris of thermal « ficients of 
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of porous materials if some degree of control is 
to be exercised over the thermal coetlicicnt of 
permeability. Since the degree of “ roug mess’ 
of the internal surface apparently determines this 
coellicient, minimization of the variation of 
permeability with temperature might be achieved 
by controlled particle sizing, more gradual sinter- 
ing, ete, 

Pore size distribution measurements of sintered 
porous media are valuable for the determination 
of the existence of structural similarity among 


various porous materials. For different sintered 


compacts which are composed of particles of the 


same shape, e.g. spheres, the mean pore diameter 


caleulated from these measurements gives an 


indication of the order of their permeabilities, 
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NOTATIONS 


cross sectional area of porous core in the 
direction of flow 

average pore diameter 

gas mass velocity (based on empty 
cross-section) (ib her ft*) 
conversion factor (32-17 ft by see*) 
permeability (maar) 
length of porous sample (in.) 
molecular weight 

Reynolds number 

constant 

pressure (p.S.1.a.) 
volume tlow rate at 20 C and 14-7 p.s.ia. 

(s.c. fh.) (based on empty cross-section) 

yas constant 

yous temperature 

vas velocity 

linear distance co-ordinate 

compressibility factor 

viscous flow coeflicient 

inertial flow coeflicient 

a difference 

absolute viscosity of the uo (Ib ft hr) 
density of the gas (Ih ft*) 
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Longitudinal mixing in packed beds 
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Abstract Longitudinal mixing properties were investigated for a water stream flowing through 

und 4in. tubes packed with 1-3, 3-0 and 3-2 mm spheres over the Revie kis number range 

to 4500 The method of injpectu step function of salt solution thir ugh «a simulated plane 

source was used, coupled with electrical conductivity measurements with very small probes 
This method is not lev the ssutnption fa that velocity profil 


The statistion ‘ ‘ by ts tested and found to ible over the 


entire range of | bles s 


tuclied lhe results are reported in terms of an eddy diffusivity as a 
function of hvdraulic radius and line velocity which are proposed as the characteristic length and 
velocity. terms 

The observed eddy diffusivities c somewhat lower than those of other investigators who 


used radially integrated rather tl int measurements of the breakthrough curve It is shown 


that the integrated results are t hivh due to a non-flat velocity profile. 


Résumé — bes te o les proprictés duo melange longitudinal dans un courant 
passant dans des tubes de 2 a 4 tn immis de spheres de 1. 5.00et 3.2 num pour un 
de Revnolds compris entre 4 methode consiste a injecter par impulsion une solution 
saline au moven d'une source plat Cette méthode n'est pas limitée par Phypothése dun profil 
plat cde vitess« 

Le model ‘ r Kinsrern est examine et les auteurs ont trouve qu'il 
sappliquait a | ‘ ites les variables étudiées. Les resultats sont donnes au moven 
dune diffusion turbulente for tu ravon hyvdraulique et de la vitesse lneaire qui sont proposes 
comme termes caracteristiques ce et de vitesse. 


Les valeurs de la diffusion rbulente observee sont un peu plus fables que celles dautres 
chercheurs. Ceux-ci utilisérent des valeurs movennes prises sur le ravon au lieu ce valeurs 
measurees par pomt sur la 


Il semble que les valeurs moyennes sont trop élevées du fait d'un profil de vitesse non linéaire 


Zusammenfassung derheiten der Laingsdurchmischung wurden fiir cinen Wasser 
strom untersucht, der durch cin m und 10 em Rohr fliesst : das Rohr ist mit 1.5; 3.0 und 
3.2 mm Kugeln gefallt fir emer tereich der Revnoldszahl von 3 bis 4 500. Die Methode des 
schrittweisen Einspritzens emer losung durch eine sinmulierte, ebene (thichenfoOrmige) Quelle 
wurde angewandt, zusammen mit elektrischen Leitfahigkeitsmessungen mit se hr kleinen Probe 
mengen. Diese Methode wird nicht dureh die Annahme eines ebenen Geschwindigkeitsprofils 
cingeschriankt 

Das von Kinstein vorgeschlag statistiche Modell wird geprift und als anwendbar fiir den 
vesamten Untersuchungsbereich ler Variablen gefunden Die Ergebnisse werden in Form 
eines turbulenten Diffusionskoeflizienten als Funktion des hyvdraulischen Radius und der linearen 
Geschwindigkeit dargestellt, cdi s charakteristische Lange und Geschwindigkeitsausdriicke 


vorges hhage n werden 


Die beobachteten turbulente Diffusionskoeflizienten sind etwas kleiner als die von anderen 
Autoren. die radial integrierte und keine Punkt-Messungen fiir die Durchbruchskurve verwandten, 
Es wird gezeigt, dass die inteyrierten Ergebnisse infolge cines unebenen Ges hwindigkeitsprofils 


zu hoch sind 
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INTRODUCTION 

Dining the last few vears, increased attention 
has been focused on the problem of longitudinal 
mixing in packed bed systems, using both gases 
and liquids as the fluid) mecha. This mixing 
phenomenon is important because of tts miluence 
on the performance of packed bed reactors, 
mixers, and extractors. In the case of reactors 
and extractors, this mixing process decreases thy 
effectiveness of the apparatus. A study of the 
extent and variation of the longitudinal mixing 
properties of pac ked bed systems leads to a mor 
thorough understanding of the (haid dynamics and 
tubulence properties of such flow systems. 

Several investigators [1, 2. 3, 4.5, 6,7, 8, 15, 24 
have studied longitudinal mixing im the liquid 
phase in packed beds, using various types and sizes 
of packing material. The most popular methods of 
tracer injection have been the pruilse function |} 2. 5 
and the sinusoidal concentration wave 
7. 8). All of the experimental results have been 
analysed on the basis of the diffusion model. In 
the case of the diffusion model, the familiar type 
of differential equation is solved, usually assum 
ing an infinitely long packed bed. Some investi 
gators [2, 6, 7, 9] have solved equations represent 
ing the longitudinal mixing process as analagous 
to that of dispersion resulting from a fluid stream 
passing through a series of stirred tanks. When 
the results of the stirred tank model are mace 
consistent with the diffusion model, it is possible 
to estimate a value for the longitudinal Peclet 
group. 

The purpose of this investigation was to 
determine the longitudinal mixing characteristics 
of a packed bed, using expe rimental methods 
which allowed maximum possible freedom from 
previous assumptions limiting the validity of th 
results, and to test the applicability of a statistical 
model. This model does not require the diffusion 
assumption which ts valid only in long beds nor 
the x rfect cell mixer assumption which contains 
parameters diflicult to evaluate; nor does it 
require the introduction of an additional para 
meter, such as mixing elliciency, or a bed capaci 
tance term. The statistical theory was first 
introduced by Eixsrrein [10] to describe the 


motion of pebbles under the influence of a water 
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stream. It has the advantage of being applicable 


for anv length of bed, however short. 


THEORY 
The theory for longitudinal mixing in packed 
beds based on the ErinsTrern statistical model has 
heen derived and discussed in detail elsewhere 
11. 15, 


discussion, the essentials are summarized. The 


For convenience in use for the present 


model considers the motion of corpuscles ot tracer 
material in the mainstream fthaid. The tracer 
corpuscles alternately experience rest phases and 
motion phases, whereby the motion phases require 
a small amount of time compared to the time 
interval required for the rest phases. The rest 
phases may be issumed to occur at stagnation 
points at each packing particle and the motion 
phases occur in the void spaces through which 
the fluid moves at high velocity. 

The probability density function for the 
corpuscle in’ time and space (one dimensional, 


longitudinal) is given by 
dt exp | r) dzdr (1) 


where p(2, 7) is the probability density for the 
tracer corpuscles at a longitudinal distance z from 
the point of introduction and time ¢ after the 
introduction, 

is the dimensionless longitudinal distance from 
the tracer source, 

r is the dimensionless time since the corpuscle 
of tracer entered the system. 

As described previously [11, 15], equation (1) Is 
applied to a consecutive set of motion phases 
and w rest phase s after which the tracer corpuscle 
at the time 


The considerations are treated in such a manner 


is found at the position 


that all possible paths for arriving at N at time 
T are taken into account, finally arriving at the 
eXpre 


Fy CXp | N T) (24 NT) dT (2) 


a dimensionless time. 


W it 


\ a dimensionless longitudinal distance 


(from the tracer source) 


12 
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e Bessel function of zero order and first 


kind of an imagmary argument. 


The symbol Fy .,, represents the probability of 


N at 


where tracer corpuscles orgimate (at 


finding a tracer corpuscle at position 


0) continuously and at a steady rate starting 


time 
at zero time. This definition of Fy corresponds 
to a step function tracer input. 

The 


ment were : 


important assumptions in the develop- 


1. Flat mainstream velocity protile over the 
small region of interest, 
2. The probability of a tracer corpuscle 


coming to rest is independent of time and position, 


3. The time required for a motion phase ts 


small compared to that required for a rest phase. 


Equation (2) is difficult to evaluate ; therefore 


an approximate form [12], [15] may be used: 


The relationships between T and N and the 
more familiar eddy diffusivity FE and time f 
are given by 
ZU 
\ - (4) 
7 (5) 
(0) 
Z 
N : d, (7) 
Z 
Equation (4) is a detinition of the eddy 
diffusivity in terms of the Einstein statistical 
equations. The numerical value of the eddy 


diffusivity is equal to that obtained from the 


same data analysed on the basis of the diffusion 


model when the beds are long enough (N 50). 


The familiar breakthrough curve for a step 


function input is obtained by plotting Fy ,,, as 
in equation (3) versus (7 N) for a given value 


of N. 


The most convenient method for evaluating 


Carens and 


J. 


M. 


the Peclet group from the F curve ts to obtain 


the slope ()F° )@) at the point @ 1, sinee from 
equation (2) one may obtain [11]: 
MF 1, (2N) 
y fo (8) 
For large N. it may be shown that 
ad d 
pe 
Z M0 Z 
Equation (8) indicates that this method of 


calculation of the Peclet group is very rapid and 
requires considerably less calculational procedures 
12, 13). 


than other methods [2, 4, 5, 8, 9, 


EX PERIMENTAI 


The method employed for investigation of the 
axial mixing properties in packed beds was that 
of a step function injection of tracer solution. 
In this the 
solution the 


conductivity, 


nitrate 
that 
maimstream 


tank 


case, tracer was a sodium 


method of 


Phe 


constant-head 


and detection was 


electrical wis 
from a about 


the 


tap water 


35 ft above test section. 


ie ral flow 


apparatus 


\ schematic flow diagram of the 
Fig. 1. 


apparatus Is 


A 5 


tise dl 


shown in hip. regenerative 
the 
of pump is capable fo deliver- 
at 


fluctuations. 


typ 


turbine was for maimstream, 
since this ty yn 
with 


Calibrated 


ing medium flow rates high pressures 


negligible surges and 


rotameters and globe or needle valves were used 
the 


was designed 


The thow distributor which 


discussed in detail elsewhere (ut). 


flow. 


to 


to regulate 


wive a velocity protile Is 
The test section for most of the data re ported 


here was a two-feet high, 2 in. glass pipe section 
composed of four 6-inch lengths tlanged together. 
the conductivity cell to emerge 

test between gasketed 
Son bin. 
a different prolx arrangement 
The tank 


allows drainage of the thud without influence on 


leads 
the 


from a 


allowing 
the 
flanges. 
test 
described 


section 


data 


from 
are included 
section and 
elsewhere overtlow 
the test section. 

The packing particles were spherical in all 
1°33, 


cases. The particle diameters studied were 
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340 and 82mm. Bed heights varied from 5-6 


to 


over flow tank 


Probe holder 


Hypoder mic 
needies 


Schematic flow diagram 


clor syste m 


The tracer myector system consisted of a 
piston and cylinder driven by a Graham variable 
speed drive and gear reducer unit. The speed of 
the piston was determined as a function of the 
drive speed regulator position. 

The injector head for the packed bed was 
constructed from 18 gauge hy pod rmic needles 
arranged on a % 82 in. square lattice. There wer 
37 needles in the Zin. plane source, and 156 
needles in the 4 in. plane source. 

Calculations were performed in order to 
determine a reasonable needle spacing for the 
assumption of a plane source. KLINKENBERG 
et al, (12) solved the equations describing the 
spreading of a tracer from a pomt source into a 
moving fluid stream where axial and radial 
mixing rates are not assumed equal, KLInKEN 
neRG’s results were used for calculation of the 
lateral concentration profiles at various distances 


downstream from the source. The method of 


superposition was used for the construction of 
concentration maps for various needle spacings 
at various distances from the source. Details of 
the calculations are described elsewhere {11}. 
The results for these injectors show that a distances 
of about yin downstream from the injectors 
there is only 9 per cent “ ripple in the concen 
tration map. This shows that the radial diffusion 
necessary to establish a plane concentration 
front occurs very rapidly with respect to the 
longitudinal diffusion which is detected at a 
distance of 18 in. downstream. Hence, the radial 
diffusion effects have a negligible effect on the 
longitudinal diffusion results. 

Fig. 2 shows a photograph of the 4 in. plane 
source, The body of the source is hollow, allowing 
the tracer solution to flow into the ends of the 
Lin. cross tubes supporting the needles. Caleula 
tions showed that the pressure drop through 
each needle was about 100 times as large as the 
pressure drop through the § in. cross tubes, and 
much larger compared to the pressure drop 
through the inside of the injector body. Therefore, 
an even distribution of injectant was to be 
expected. This was obs rved to be the case. The 
injector was mounted directly above the flow 
development section as shown in Fig. 1. 

The solenoid valve arrangement allowed a very 
rapid action either for turning on or turning off 
the flow of the injectant into the test section. 
Two valve arrangements were used to test the 
effectiveness of the step function generation by 
the valves 

The first unit was constructed such that there 
was a General Controls Maenetic Valve typ 
K21-F within tin. of the injector head, and 
another on the by-pass pipe as shown in Fig. 1. 
By throwing a single switch, one valve was 
ope ned and the other was simultaneously closed. 
The valves were installed such that the direction 
of tluid flow assisted fhe rapid closing action. 

The second unit was composed of a General 
Controls type K21-F Magnetic Valve on the 
by-pass pipe, and a Jamesbury stainless steel 
rapid-action ball-type valve actuated by two 
large solenoids and an aluminium linkage system 
arranged to give a push-pull action for turning 


the valve stem which had a long aluminium 


~ 
" ‘ 
| | 
onductivity 
probes 
Test sectior 
ect A 
. ect ond 
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< 
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valve was 
head. 
the 
The initial 
made use of a 
Ihe cell 


connected through the electronics to a high speed 


lever mounted on it. The Jamesbury 
in. 


mounted about from the imyector 
Both 
Peclet group at identical flow rates 
of the 


cell at a 


valve units identical results for 


testing valve units con 


ductivity needle WAS 


tip 


recorder having chart "5 mM sec, 


Both valve 


function response even down to very 


speed afl 
arrangements gave a very sharp step 
thon 


rates (about 0-0006 eal min) 


Probes and clectronics 


for most of the 
Fig. 3. 
structed from packing units of a size 
to the 


were 


The probs used pote ke ad bed 


studies is shown in This probe is con 
identical 
Ihe 
that 


‘ j weet t he 


rest of those used in the bed leads 


strung with perforated beads so they 


could bn led out of the flanges 


with a minimum of disturbance to the rest of the 
The 
parallel and the plane of the probe w 
dicular to the 


that 


electrode faces of the probe were 


packing. 


is 


direction of nd flow im 


within the 


order no hold-up be allowed 


Packed bed conductivity cell 


and 


J. M 
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cell, The 


manner that 


only the 


constructed im 
in the 


probe contributed 


conductivity was 


such a inter 
stitial 


to the 


Space formed by the 
conductivity. 
constant 


used as the 


The basic circuit is shown in Fig. 4. 
20 ke 


was connected across the 


voltage oscillator voltage 
conductivity cell 
The 
across ris proportional to the conductance of the 
eell, 1 OR, 


The relore 


signal il 


soures 


and the current detector resistance } voltage 


solution in the providing r is always 


small compared to R. by amplifying 


and demodulating the voltage 
directly proportional to the solution conductance 
reco} cle al. SIN 


that the 


may be parallel channels were 


‘ ells 
The amplifier 


conductance of up to 


SIN could bye recorded simultaneously 


demodulator channels had a linear 
The 


varnibl 


effects tor 


vain of up to vain Was con 


tinuously from 0 to Capacitance 
both types ol cells were measured and found to 
Careful 


utilized 


be neghoible. shielding vroundine 


trees were 


hu t. Block diagram for 


shows a block diagram of the overall 


Fig. 4 


eireuit for one channel. Six channels in parallel 


were used for most of the work. The pip box sa 


electronically 


trace of the 


device which placed a pip in the 


reference recorder every 1 60th of a 


second, The box was activated by a 20 ke signal so 
that the 
The pip box also puts a separate additional pap 


chart at tha 


error in pip placement is 
when the 
The 
start pip placement was within 


to the 


on the moment myector 


solenoid valve switch was thrown 


of the 


This slight error is duc 


small but variabk 


delay in the solenoid action due to the phasing ol 


the G0¢ voltage source for the solenoids at the 


moment that the switeh was thrown. This error 


q 
VOI 
12 
constant 
voltage source 
Ampiitier 
7 lator 
nsuictec Sac OF er 
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Switch mo ~Pip 
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was negligible except at extremely high Reynolds 
numbers, 

The high speed multichannel recorder was a 
Minneapolis Honeywell Visicorder model 906, 
The traces on the chart were the result of the 
motion of mercury vapour light spots reflected 
from the mirrors of the bank of galvanometers 
mounted in the instrument. The deflection of 
the light spots was very nearly proportional to 
the voltage across the galvanometer and damping 
resistor. The frequeney response of the instru 
ment is such that transients up to 300 ¢ s may be 
followed with an accuracy of about 2 per cent. 

Detailed circuit diagrams and descriptions for 
all of the electronic equipment except power 


supplies ure available else where 11 


Procedure 


The injector drive setting was determined 
such that the linear speed of the injectant leaving 
the needles was equal to the average linear speed 
of the mainstream as calculated from the tlow 
rate and fraction voids. 

When steady state injection had been attained 
the run was started by turning on the chart 
drive and then causing the mnjectant stream to bn 
instantaneously diverted to the bypass tube. 
At the same instant the pip box automatically 
placed a sharp start pip in the reference trace on 
the recorder chart. 

The response curves for the shutting off of the 
tracer were used because of the advantage in 
being able. to pre-set the gains and knowing 
the relative positions of the traces. It was found 
to be easier, in general, to turn off the tlow of the 
injectant stream very suddenly than to turn it on, 
due to volumetric capacitance and momentum 
effects in the injector body. With proper care, 
however. identical results could be obtained from 
both types of curve. Data wer collected for 
low rates corresponding to Reynolds numbers 


Ud. p en from about 3 to about 4500, 


REDUCTION or Data 


In order to develop a rapid, simple method of 


calculation, an exact expression having the form 
ofl approximate equation (9) is desire d. However. 


it is not possible to determine directly the location 
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of the port (T N) 1 from the data without 


resorting to trial and error. For this reason, 


the slope of the F curve ts taken at a pom which 


may be located directly and simply. Such a 
point is F 1/2. 
Knowing the slope at the pomt F 1 2. the 


slope at the pot 7T/R I be caleulated. 
using the derivative of equation (2) with respect 
to (7. N). Assuming a value for N. the corres 
ponding 7 for F 1 2 mav be calculated from 
equation (3), and these values of 7 and N are 
inserted into the derivative of equation (2 
Then at F 1 2iscaleulated. The value 
for the slope at T/R 


(8) for this same assumed N. The ratio of the 


1 is found from equation 


two slopes is detined as the correction fac tor 


+ 10) 

Therefore ()F 0 , may be obtained directly 
from (JF d@)p_,. from a_ plot of versus 


(OF 

Since equation (9%) ts only approximate, it is 
desirable to have a relationship between an 
equation having the form of equation (9), and 
the exact solution. The relationship ts ¢ xpressed 


in the form 


\ | 4) a 
| exp (2N 


For a general detinition of x. see the App nadix, 
Equation (11) defines x for use in the present 
calculations. The calculation of the Peclet group 
now follows the form of equation (9). but is exact, 


whereas equation IS approximate 


Z 

Since the electrical conductance of the sodium 
nitrate solution which was used as a tracer ts not 
exactly proportional to the concentration of the 
sodium nitrate, additional correction ts necessary 
in order to prov ide a means of direct calculation 
of the Peclet group from the slope at F 1 2on 


the recorder chart. The form of the empirical 


e 
25 


J. Carns and 


equation relating the salt concentration C to the 


conductance L is: 


( al! (13) 
where a constant (strong function of tempera- 
ture). 
h constant (weak function of tempera- 
ture). 
Kaquation (13) was used over the 
4 ( 2 N, with good precision. 


The correction factor (which is developed in the 


Appendix) takes the form : 


which ts a weak function of temperature. 


It was found that even though (1 Ls) was 
verv nearly equal to (y y) (where y is. the 
displacement of the trace on the recorder chart, 
and the subseript indicates the steady state 


value during injection of tracer material), 


| (15) 


L L } 12 


Introducing equation (14) and equation (15) to 


equation (12), the tinal working equation for the 


determination of the Peclet group becomes 


ly pron 
Z da\ \* 
\dtl 
d yu 
where slope of the breakthrough curve 
da 
on the recorder chart taken at 
1/2. 
L steady state conductance valu 


during injection. 


line velocity of the 


U, dz where 


pont velocity 
dy dt chart speed, 
A? — pt a2, 


The term in the large brackets is (dF dQ), , 


J. M. Prausnrrz 


Equation (16) provides a rapid method for 
reducing data from the experimental step 
function method. For a plot of A? 1 versus 
(XF 0); SOC Fig. 5. 


discussion are given in the Appendix. 


Details of the foregoing 


Pia. 5. Correction factor for obt ining De 
The eddy diffusivity E is found from 
Ud UZ 
/ (17) 

Pe 


where U is the characteristic fluid velocity. 
Fig. 6 was prepared as a test of the validity 
of the model. The eX rimental points were 
obtained from the recorder chart. The linear 
velocity in the term U¢% Z was found from the 
start pip on the chart, the known Z, and the time 
required to reach the concentration which 
corre sponded tor the point for (7 /N) 1 from 
equation (3). This is a trial and error procedure 
lor (T Vv) 1 is very close to 


1 2. so only one trial was necessary for suflicient 


however 
accuracy percent). The agreement between 
expermment and theory is well within experi 
mental error. 

In reducing the data as outlined above, the 
question of characteristic velocity arises. Previous 
investigators |2. 7, 15, 16, 17) have used the inter 
stitial velocity as calculated from flow rates and 
fraction voids U, « This is only satis 
factory when the assumption of a flat velocity 
profile is valid [11, 18}. This assumption is very 
seldom justitied, as may be seen from results on 
velocity protiles in packed beds | 1S, 19) and from 


a discussion which follows in the next section. 


| 
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Comparison of experiment to theory 


Curve: F 


Data points: run 5, set 3, cell 13, Re 


The characteristic velocity which is thought to 
be the most valid here is that obtained from the 
the tracer through the column 
as characterized by the 

the start pip and the 7) N 
closely corresponds to the F 


rate of motion of 


time interval between 
1 point, which very 

1 2 point in the 
range of variables studied here. This velocity ts 


called the line velocity and is defined as: 


the distance from the 


injector to the point of observation. 


longitudinal 


where Z 


the point velocity. 


The nature of this definition is such that the 
velocity is averaged along a longitudinal element 
of the but 
velocity variations have little or no influence on 
the 
This detinition is also advantageous since a point 


tube, not radially, so that radial 


line velocity as measured experimentally, 


velocity is not used. The velocity is averaged 


over essentially the same path as that corres- 


sy T sy N 
128 Z 


24-0 in. in. sec 


ponding to the longitudinal mixing process, and 
the 


by this velocity than by any other. 


hence mixing is characterized more closely 


ResuLTs 


In order to summarize the results it is necessary 
the affect the 


parameter. Previous investigators 


to know how variables eddy 


diffusion 


7, 15) have plotted the Peclet group 


superticial velocity. 

fraction voids, 

particle diameter. 

kinematic viscosity. 

eddy diffusivity as calculated from 
the 


using 


experimental results 


the 


various 


€ as characteristic 


velocity. 


0-6 
0-2 
JOL. ( 
] 9 ae 0:90 )-95 00 10 20 
O40 
| 
2, 5, 
yp Pe od, 
versus the Reynolds number 
where 
d 
4 
E 
a 
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One investigator [2] plotted FE versus Re, and 
The effect 


of viseosity on eddy diffusivity has been inves- 


another plotted versus Re. 


tigated in the range of from 0-95 to about 28 cp 
by Eeacn and Wurre [5]. 
table effect. 


have been proposed predict any typ of viscosity 


There was no detec- 


None of the models or theories which 


dependence, 

The model proposed by PRAUSNITZ [233] cIVes 
a proportionality between the longitudinal eddy 
diffusivity and the product of interstitial velocity 
The 


in this work also predicts a 


and particle diameter. EINSTEIN statistical 


roc | propose d 


proportionality between the longitudinal eddy 


diffusivity and the produet of the characteristic 
velocity (UC) proposed in this work) and a charac 
teristic length term. Therefore it does not appear 
eddy diffusivity versus 


appropriate to plot 


Revnolds number since the eddy diffusivity is 


hho way cle rk nt on the kinematic viscosity. 
The 


proposed here as 


hvdraulie radius as extended by 
a characteristic length 
because it is not 


term for plotting 


to parcke al 


Purpose 


limited beds, whereas the particle 


PRAUSNIT? 


diameter as a characteristic length term loses its 
significance when applied to anything except a 
packed bed. 

For data summarizing purposes, it is proposed 
that FE be plotted against inl’, «. where m is the 


hydraulic radius term given by 


volume of thiid (free volume) 


wetted area 


In the limit of the open tube, 4m ts equal to the 


tube diameter. For spherical particles im a 


eviindrical tube, it can be shown that 


De 


"3 2(D 4) (1 


where D tube diameter. 
é fraction voids. 


particle diameter. 


For the 32mm spheres in a Zin. tube with 

O-3T9, 4m 
Fig. 7 shows the results of this work compared 

with the The 


ordinates for Fig. 7 that the 


1-217 mm. 


results of previous workers. 


were chose Stik h 


4M UA), 


Kaddy diffusivity in packed beds (based on average 


em’ /sec 


velocity) 
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Jacques and Vermeulen 
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Kddy diffusivity in packed beds (based on line velocity ) 
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D Zin 


3-196 mum tin 


mam D bin 


2-06 mom 


independent variables can be evaluated easily, 
not because they are believed to be the best for 
characterizing the mixing in packed beds. A 
knowledge of the velocity protil is not necessary 
in order to apply the results shown. The co- 
ordinates proposed as a result of the present 
work are versus where is calculated 
using the same equations as for FE, but using if 
rather than U, « as the characteristic velocity. 
The results of this work are shown on these 
co-ordinates in Fig. 8. 

The more familiar curve using Pe as_ the 
ordinate is shown in Fig. 9 The abscissa, how- 
ever. is not Re. since it has been shown that Pe 
does not depend on kinematic viscosity, The 
term 4m’, « has been chosen for the abscissa, 
Included on the plot for comparison purposes ts 
Pe’. the 
used os the characteristic velocity, versus tml’). 


prime designating the fact that U') was 


Fig. 9 indicates that the present results lie 


D tin 


above those of other workers. This general 
discrepancy may be explained in terms of the 
assumptions made by the various investigators, 
coupled with their experimental methods. 

Some workers [2, 5, 3] used a dye tracer 
technique, coupled with a light transmittance 
technique whereby the light is transmitted across 
the whole tube cross section at a given Z level, 
which is devoid of packing. The pulse function of 
dye which was used in some work [2, 5] was injected 
into an unpacked section where flow not charac- 
teristic of the packed section prevailed. Entrance 
effects and non-flat velocity profile effects could 
Another 


pulse function method [5] was used in which the 


be important in such an apparatus. 


pulse was injected into the mainstream before it 
entered the expansion section leading to the 
packed test section. Here, entrance effects as 
well as non-flat) velocity profile effects can 


influence the results. Even neglecting the 
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entrance effects. the non-tlat velocity profile 1 method, since the breakthrough curves are 
the test sections can influence the results to a averaged across the tube before the Peclet number 
previously unsuspected extent. is calculated. Fig. 10 shows that the curve of 
\n example using a step function input best the average of the two components has a slope 
serves to explain the inaccuracies encountered at the F 1 2 point which is less than the slope 
when the breakthrough curve is determined of either of the components. This average curve 
by averaging across the test section as is done would by necessity vi ld a Peclet yroup which is 
with the light transmittance  te« hnique, For misleadingly low. The Peeclet group calculated 
simplicity, it is assumed that half of the fluid is from the resultant curve is 0-55. using U, U,/¢, 
travelling at OF tly «, and the other half resulting in an error of 45 percent. The apparent 
at U, O-o 0, «. giving an average ve locity of mixing is * better ” than the actual mixing due to 
«. The thuid with the higher speed is near the the non-flat ve locity profile. Even though the 
wall as an annulus; the slower fluid forms a example given was for a step function injection 
central core. Two observation points at some the same error would hold true for the pulse function 
distance Z from the plane source are located type of injection. Experiments such as the ones 
such that one is near the centre, the other near reported here are not sensitive to the flatness of 
the wall. Fig. 10 shows the two calculated — the Velocity protile. 
breakthrough curves for a Pe of unity. Equation Another set of investigators [15] assumed that 
(3) was used for these calculations, with N 100. a plane source was adequately approximated by 
The breakthrough curve closest to the origin seven point sources in a column which was 8 
corresponds to U, lil, «. The dotted curve particle diameters in diameter. Calculations 


is the average of the two curves which in essence — similar to those described for the Injection system 


is the method used in the light transmittance used in this work show that a distance of 5 particle 
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Fig. 10. Effect of non-flat velocity profile on averaged breakthrough curve. 
Slope of resultant (dotted) curve 209 re O55 Actual Pe 1-0 (45 per cent error) 


diameters from the source was necessary before 
the * ripple “in the concentration contour map 
was 11 per cent taken along the diameter which 
contained the least ripple. The column was 
effectively only 25-35 particle diameters long, 
and hence about 20-25 per cent of the length of 
the column was subject to conditions which were 
not described by the boundary conditions assumed, 
Such a svstem would vield low results for the long- 
itudinal Peclet number since the radial mixing 
from a point source creates apparent axial mixing, 
even in the absence of true longitudinal mixing. 
Hence. the observed resultant longitudinal mixing 
is higher than if a true plane source had been 
used, 

From the above discussion it can be seen that 
the Peclet groups of all of the investigators cited 
would tend to be lower than the results for the 
Peclet group based on a true plane source and on 
point measurements of the break through curve. 
Correspondingly, the results for the eddy diffusi- 
vity as reported by other investigators are higher 
than those reported here. 

By reporting the eddy diffusivity as a function 
of the product of the hydraulic radius and the 
average velocity as in Fig. 7. the results of this 


investigation are somewhat high due to the use 


31 


of the average velocity in the calculations ; 
but they are not as high as those reported by 
others due to the fact that the break-through curve 
was not averaged across the tube as discussed 
above, but was observed at various potnts in the 
tube. 

The presence of a non-flat velocity protile would 
be expected to have a negligible effect on the 
present results due to the fact that near the 
source, all of the injectant across the tube is 
travelling at a velocity equal to the average 
velocity. As the velocity of the tracer accommo- 
dates itself to the main-stream velocity profile. 
longitudinal diffusion is very rapidly taking 
place and reducing the longitudinal concentra- 
tion gradient. At the same time, a radial con- 
centration gradient which is expected to be small 
compared to the longitudinal gradient is being 
established. The radial eddy diffusivity is only 
110th to 1 Sth as large as the longitudinal eddy 
diffusivity, and hence mass transfer in the radial 
direction is much less than that in the axial 
direction. 

Fig. 7 also shows that the results of McHenry 
and Witue tm [7] for gaseous systems and the 
results of others for liquid systems follow one 


line whereas there is a wide divergence when 
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Sine 


One merely measures (1 2)"/° of the distance between the 


two steady state values (before and after the step function) 
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Some aspects of liquid flow through fan spray nozzles 


N. Domerowskt, D. Hasson*® and D. E. Warpt 


Department of Chemical Engineering, Imperial College, London 


(Received 5 August 1950) 
Abstract—The flow pattern of the liquid sheet produced from a rectangular-orifice fan spray 
nozzle is investigated with particular emphasis on the manner in which its thickness varies from 
the orifice to the point of breakdown. The trajectory of its boundaries is analysed on the basis 
of a simple theory. 

The results of the study can be summarized as follows : 

(a) The streamlines of the spray sheet are straight and unaffected by the curved boundary. 
The liquid velocity along a streamline is constant, everywhere in the sheet and independent of its 
Viscosity 

(b) The thickness of the sheet is at any point inversely proportional to its distance from 


the orifice, and can be expressed, for a given set of operating conditions, by a thickness parameter. 


(c) At relatively low injection pressures, the thickness parameter for a given nozzle and 
liquid of low viscosity is found by experiment to be a function of surface tension and a function 
of (Pp) 2 » where P injection pressure, p density, 7 Viscosity At the high injection 


pressures the thickness parameter ts only a function of (Pp) . 


(dy) The trajectory of the sheet is a function of the Injection pressure, sheet thickness and 


surface tension, and independent of liquid density 


Résumé—Le diagramme d'écoulement d'une nappe liquide issue dun gicleur a orifice rectan 
gulaire est examiné en mettant spécialement accent sur la variation de son épaisseur depuis 
orifice jusqu’au point de contraction. La trajectoire de ses limites est analysée sur la base d'une 
théorie simple. 


Les résultats de cette étude peuvent etre resumee aittst 


(a) Les lignes d’écoulement de la nappe pulvérisée sont droites et ne sont pas affectées 
par la limite incurvée. La vitesse du liquide le long dune ligne d’écoulement est constante partout 
dans la cappe et est indépendante de sa viscosite 

(b) L’épaisseur de la nappe est a un point queleonque inversement proportionnelle a sa 
distance de Vorifice, et peut-étre exprimée, pour une série donnée de conditions de fonctionnement, 
par un parameétre quelconque, 

(c) L'expérience a montré que pour des pression dinjection relativement basses, la para- 
métre d'épaisseur pour une tuyére donnée et un liquide de faible viscosite é est une fonction de 
la tension superficielle ct une fonction de ou P Pression d’injection, p densité, 
viseosité. A de hautes pressions d’injections, la paramétre d’épaisseur est seulement une 
fonction de 

(d) La trajectoire de la nappe est une fonction de la pression d’injection, de Pépaisseur de 
la nappe et de la tension superticielle, elle est indépendante de la densité du liquide. 


Zusammenfassung——Das Fliessverhalten der Fliissigkeitslamelle, die durch die rechtechikge 
Offnuny einer Geblisespriihdiise entsteht, wird untersucht. Besonderer Wert wird auf die 
Feststellung der Dickeniinderung von der Offmung bis zum Auflésungspunkt gelegt. Der Bahn- 
verlauf der Lamellengrenzen wird auf Grund einer einfachen Theorie analysiert. Die Untersuch- 
ungsergebnisse kOnnen wie folgt zusammengefasst werden : 


*Now at the Technion, Haifa, Israel. 


+Now at British Oxygen Engineering Co, Ltd., Edmonton, London. 
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konstant 
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Die Stromlinien der Flissigkeitslamelle sind gerade und werden durch die gekriinonte 


entlang einer StrOmiinie rst 


fiberall in der Lamelle und unabhiingig von ihrer Viskositat 


Die Dicke der Lamelle ist an jedem Punkt umgekehrt proportional zur lentfernung 


von der Offnung und kann fiir cine gegebene Zahl von Betriebsbedingungen durch einen Dichen 


Parameter dargestellt: werden 
tc) 


fiir eine bestimmte Diise und en 


spannung und eine Funktion von (P, 1/2 » ist 


Bei hohen Spritzdrucken ist der 


(da) 


und der Oberflachenspannung 


INTRODUCTION 
THE PERFORMANCE of spray nozzles has received 
the last de« 
in combustion |1 
3 


intensive study during vle because 


chemical 
In 


applications, the function of the nozzle 1s to 


of their increasing us« 


2) and agricultural engineering. most 


disintegrate a liquid stream m order to create a 


large surface per unit volume, and to ist 
with control the resulting drops uniformly im 
space. In order to obtain a greater understanding 


of the processes involved. work has been directed 


towards analysing the hydrodynamics of 
through a wide variety of atomizers, establishing 
the 


determining the 


basic mechanisms of drop formation, and 
resulting drop-size distribution. 


In recent vears, particular attention has been 


paid to the atomization of flat spray sheets 
[4-12] because of their simplicity and because 
they prov ide a convenient model for both ex 
pn rimental and theoretical study. These sheets 
can be produced in a number of ways by the 
impingement of liquid streams or by their impact 
on a_ solid stationary or rotating surtace. \ 


simple method is to use the fan nozzle, a device 
which enables the properties and dimensions of 
the sheet to be easily controlled by varying the 
approach passages to the nozzle orifice. 

While the 
and disintegration have received increasing atten- 
littl the 


which the sheet is formed at the orttice or about 


factors governing sheet instability 


tion, very s known about manner in 
its phy sical dimensions. 

The present paper describes the prin ipal flow 
characteristics of fan spray sheets prior to the 
of 


ferometry is employed to record the variation of 


onset instability and disintegration. Inter- 


Bei relativ niedrigen Spritzdrucken zeigt das Experiment 


Dicken-Parameter nur eine Funktion von (Pp ! 
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dass der Dicken-Parameter 


ve Flissigkeit niedriger Viskositaét eine Funktion der Obertlachen 


P Spritzdrock Dichte, V iskositiit 
2), 


Der Bahnverlauf der Lamelle ist eine Funktion des Spritzdruckes, der Lamell rhieke 


er ist unabhingig von der Flissigkeitsdichte. 


thickness throughout the sheet, and a simple 


analysis, is presented to account for the trajectory 
of the free edges of the sheet. 

Tur Mopr or Formation or A FAN SprRAY 
The principal methods of forming thin sheets 


of atomizers have been discussed 


na wick range 
by Downrowsk: and Fraser [6]. In the fan 
spray nozzle, two streams of liquid are made to 


nopinge behind an orifices by specially de signed 
approach passages and a sheet is formed in a plane 
rpendicular to the of the streams. The 
principal ts illustrated in Fig. 1 (a) which shows 


liquid flowing through a rectangular orifice formed 


at the end of a rectangular tube. Under these 
conditions the tlow through the orifice is con 
stricted in only on plane and the streamlines 
converge to form a region of pressure behind 


the 
freely spreads through the orifice limited only 
by the 
sheet can further be increased by extending thi 
to the of the (Fig. Ib). 


1 (c) shows a commercial nozzle, designed on 


orifice. A flat sheet ts produced as the licguid 


side walls. The spreading angle of th 


opening sides orice, 


Fig. 
this principle, which has been used in the present 
work. It is made of porcelain and contains a 
rectangular orifice which is produced by the 
interpenetration of two rectangular slots. 

In the absence of surface tension, the edges of 
the sheet would travel in straight lines from the 
orifice so that a sector of a circle would be formed, 
However, as a result of surface tension, the edyves 
contract and a curved boundary Is produced us 
\ typical 
flash photograph of a low viscosity fan spray ts 
It the 


the sheet deve lops bevond the orifice. 


illustrated in Fig. 2 (a). shows that as 


| 


Sheet 
Region of 
hgh pressure 
—_ 
(a 
7 iquid 
heet 
4 
Region of 
hgh pressure 
; 
| 
liquid in each edge moves along the curved 
boundary, the latter becomes disturbed and 


manner somewhat 
When 


resulting drops sustain the direction of tlow of the 


disintegrates in a analogous 


to a jet of [17]. this occurs. the 
edge at the pot at which the drops are formed, 
and remain attached to the receding surface by 
which rapidly disintegrate into 


thin threads 


streams of drops. Corresponding further to the 
characteristics of jets, the breakdown of the edges 
is restrained by high viscosity. At higher injection 
velocities (Fig. 2b), the contraction is less pro- 
nounced, and the placid sheet eventually becomes 
ruffled and disintegrates before the two edges 


malesce. 
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* Orifice orea=0-0784*0 0343cm 


Mechanism of flow through fan spray nozzles. 


Analysis of flow in sheet 


streamlines, photographs have been taken of spray 


In order to examine the nature of the 


sheets containing suspended aluminium particles. 
Fig. 3 
liquid within the sheet and as it passes through 


demonstrates the direction of flow of 


the edges into the ambient atmosphere. It was 
obtained by taking a photograph with a flash 
each 
Three 


points of interest can be observed in this photo- 


of approximately I msec duration so that 


aluminium particle traces its own path. 
graph. First, the sheet streamlines are straight 
and unaffected by the curved boundary. Secondly, 
the drops leave the edges tangentially at an angle 
different from that of the streak, and thirdly, 


[ 
4 
) 
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many of the drops at the base of the sheet cross 
over from one side to another. 

The velocity distribution along each streamline 
has been determined by taking two mu rosecond- 
flash photographs separated by a short and known 
The 


thus 


interval of time on the same plat 18}. 


movement of each alumimium particle 
provides a measure of the local liquid velocity. 
Measurements from photographs, su h as Fig. 4, 
of various sprays injected at atmosphe ric and 
sub-atmospheric air densities, have shown that 
the distance travelled on the sheet by individual 
particles in unit time is constant, its value being 
independent of either viscosity (up to 
or of the ambient 


atmosphere, Comparison 


between the velocity of the sheet as represented 
by the tracer parti les and as measured dire tly 
by the movement of the centres of perforations 
appearing on a sheet sprayed in vacuo (Fig. 4), 
or when soluble oil is added to water [6), showed 
Thus, it 


that the 


no measurable difference. may be 


concluded from these results stream 
velocity is constant along the sheet, its absolute 
value depending only on the differ ntial injection 
pressure. 

It has been previously assumed [19, 20, 21) 
that when a spray sheet develops, the increase in 
occurs at the expense of ts 


surface energy 


kinetic energy. Expe rimental evidence suggest- 
ing this magnitude of ve locity reduction has been 
obtained by WEIN BERG [22] for the con al sheets 
produced from swirl spray nozzles. His tindings 
cannot however be accepted without reservation 
that the 


the total head tub employed for measuring local 


since his results indicate diameter of 
velocities was larger than the thickness of the 
shects at the leading edge and thus was not fully 
immersed in the liquid. 

In order to clarify this point further, a high 
speed cine-film was made of the liquid sheet 
issuing from a swirl spray nozzle. Measurements 
made of the movement of local irregularities and 
waves confirm the above conclusions that con- 
trary to previous assumptions, the stream velocity 
is in fact constant throughout the length of the 
sheet. 

This finding is easily explained since no external 
force acts on the sheet and the flow through each 


N. D. Hasson and D. EF. 


Warp 


liquid streamline must satisfy the equation of 


conservation of momentum. 


this result to the equation of conservation of 


energy further predicts that the energy absorbed 
mn producing the larger surface area is derived 
from the enthalpy of the liquid. However, this 
cannot be easily confirmed since the predicted 
temperature drop along the sheet, calculated on 
this basis is too small for accurate measurement. 
The general form of the energy equation of a 
flowing fluid may be formulated as [23}. 

dq dV du d(i1,2mV*) (1) 


d 


where q net heat absorbed from the surround- 
ings. 


iW net work done by the fluid on the 


surroundings 


mg: — potential energy of mass m of fluid 
at height z above datum level. 
1 2mV? kinetic energy of the fluid, 


i total enthalpy of the liquid of mass 
m which is here detined as including 


the total surface energy. 


For a liquid sheet ejected from a orifice the 
equation may be simplified as follows, Since there 
is no change of velocity and assuming that no 
transferred to the 

0 and dW 0. 


can be shown that the potential energy change 


energy 1s atmosphere, 


d(1,2mV*) Furthermore, it 
is negligible compared to the surface energy 
that d (mgz) Oo. If it is 
assumed that the tlow is adiabatic, dg = 0. 


change so further 
Hence 

dH 0 
Since the flow is incompressible at constant 
pressure 


di =dU =0 (2 


where U is the total internal energy of the liquid 
mass m, defined as including the total surface 
energy of the liquid. 

Now U may be expressed in the following form 
(24) 


U=E 4 (3) 


The application of 


VOL 


12 


Avads uny 


3 
\ 0 4 . bd . 
A 
> 
Sey 
~* 


Fig. 3. Streamlines in fan spray sheet of water. (Exposure time, 1 msec 
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Fic. 4. Double-flash photograph of 1 per cent aluminium particle suspension sprayed in vacuo. 
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Where E — internal energy of the liquid when 
its surface area A 0. 
the free surface energy per unit area, 
or surface tension per unit length, 
T — absolute temperature, 
and the expression {y -T dy is 
dT 
the total surface energy of the liquid 


per unit area, 


Over moderate temperature changes, the surface 
tension decreases linearly with increasing tempera- 


ture according to the relation 


a LT 
where a and b are constants, 
Hence 
a constant 
d7 


Thus from equations (2) and (3), 
dE adA 


The change in internal energy can also be 


expressed as 


dE mcedT 


where ¢ is the specific heat of the liquuid, 


Therefore 
medT 
i.e. the temperature drop along the sheet is given 
by 
2A 
AT =~ 
me 


since the liquid mass issuing from the orifice has 
a negligible surface area. 
Now the increase of liquid surface per unit 
mass along any element is given by 
A 2 
m ps 


where S is the thickness of the sheet and p is 
the liquid density. 
Therefore 
2a 
AT (4) 
S 
From measurements of sheet thickness made 
by Srraus [4], a particular nozzle spraying water 
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at 20 p.s.i.g. produces a sheet which is 5-8 10-4 
em thick at a distance of 3 em from the orifice. 
At room temperature the value of the constant 
‘a’ for water is equal to 115-88 [24]. Equation 
(4) predicts a reduction of sheet temperature of 

Although the converging streamlines in the 
nozzle cause a region of high pressure to be formed 
behind the orifice, it will be assumed for the 
following analysis that the liquid flows from the 
nozzle as if there were a line source of high pressure 
perpendicular to the sheet. Furthermore. it will 
be assumed that the contraction of the edges by 
surface tension does not affect the flow pattern 
of the sheet, ice. the liqguid corresponding to the 
* vanished ” part of the sheet lving between the 
theoretical linear and actual curved boundaries is 
concentrated at the curved boundary. \ diagram 
of this simplified flow pattern is shown in Fig. 5. 

Consider an element of the sheet. of length dr, 
subtending an angle d8 at the orifice, and at a 


distance x from the centre of pressure P, 


If V liquid velocity, 
and § the thickness of the sheet. 


then dQ, the volume of liquid flowing through the 
element df in unit time is given by 


dQ = 2.dB.8.V. 


dQ 
dBaV 
Since the liquid streamlines are straight, and the 
stream velocity is constant throughout the sheet, 
dQ/dB8 and V are constant, and thus 

Ss K (5) 

where K is a dimensional constant which uniquely 
represents the variation of sheet thickness with 
distance from the nozzle. 

Because the sheet edge is curved by surface 
tension, fan spray sheets cannot be characterized 
by a spray angle. However, 20,7, the sheet 
angle at the orifice, before it is affected by the 
exit passage or surface tension, can be derived 
in terms of the thickness parameter K. 

The actual area of flow is represented by the 
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Traejctory of the sheet edge 

The spray pattern of a fan spray depends, 
among other things on two factors. the disinte- 
gration of the leading segment of the sheet 
(Fig. 2b), and the disintegration of the edges of 
the sheet (Fig. 2a). The former is predominant 
at most injection pressures used in practice 
where a high relative velocity exists between 
the liquid and ambient atmosphere. The latter 
becomes important at low relative velocities. 
Furthermore, it will be shown in a later section 
that at low intection pressures the tlow of liquuid 
in the sheet appears to be influenced by the 
contraction of the edges as a result of surface 


tension. A knowledge of the factors influencing 


the trajectory of the sheet edges should help 


towards an understanding of both the mechanism 
of disintegration and the tlow changes to which 
the contraction might give rise. 

An approximate expression for the trajectory can 


be derived for inviscid flow in the following manner. 


i 
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Trajectory of sheet edge 


In Fie. 7. Ois the origin of the sheet, @ is the 
angle at which the sheet edges first issue from the 
orifice, « the distance of a pomt 4 on the edge 
from O, and a dx, its position after a small 
interval of time dt. The chord AB is taken as 
the tangent at A since dz is very small. Now let 
a particle of liquid move along the line OA with 
radial ve locity 1. This velocity can be resolved 


into two components, one Vy, tangential to the 


sheet edge at point A, at an angle 6 to V, and 
the other normal to 
Thus 
sin 3 (7) 


When the particle reaches point A, it changes 


direction along the tangent. The rate of change 


of momentum per unit length is equal 


to the surface tension force at the free edge 
2 y per unit length. 
Le. 2y p 
or since 
S A x (equation 5) 
2ya 
= (5 
pk 


The radial velocity V. can be expressed as 
(vy vgP po where P is the differential injection 


pressure and C, is the coeflicient of velocity. 
Measurements made from double-flash photo- 
graphs of aluminium particle suspensions show 
that the 
approximately equal to its discharge coeltlicient. 


the velocity coetlicient of orifice 1s 
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Thus, the area coeflicient is equal to unity, and 
equation (8) can be rewritten as 


sin 3 (9) 
pk 


Now 


sin 
[ (da)? 
which on substituting in equation (9) and integrat- 
ing gives 


sin (8 + @)] (10) 


for boundary conditions where 8 7/2 0, 


when — 0, 

From this simple analysis of the tlow it can be 
that Or. the sheet the 
the trajectory of the 
predicted from a knowledge of the sheet thickness 


seen angle at orice, 


and sheet may be 


as expressed in terms of A. 


Meruop or MEASURING THE THICKNESS OF THIN 


Fitms or Liquip 


Thieknesses of spray sheets have been measured 


previously by Straus [4] and DomBrowsk! and Fraser [6 


over a limited range of conditions. They measured the 


volume of liquid passing through a small section of the 
calculated the that 
the sheet velocity could be expressed by Cy P 


sheet) and thickness by assuming 
This is a simple method, but suffers from the disadvantages 
that the the 
the sheet, that the liquid velocity must be assumed, and 
that the 


representative 


insertion of flow measuring device disturbs 


direction of flow must be known for taking a 


sumple. A number of indirect methods 


have also been employed. 
Dixon et al. [12] and Srraus [4] measured the variation 
of light absorption through dyed liquid sheets. DuKLEr 


and BrrorLin 7 employed an electrical capacitance 


method for measuring the thickness of thin liquid films 
[14 


thickness of 


flowing down glass tubes has patented a 
the 


materials using the absorption 


method of mensuring thin sheets of 


of B-rays. These methods 
suffer from their inability to detect the small changes of 
thickness 


require complex equipment and they necessitate a large 


which occur in a fan spray sheet, they may 


number of measurements to obtain the thickness over the 


entire sheet 


Since it has been shown [6] that smooth sheets can be 


obtained from spray nozzles of the type shown in Fig. 1 (c), 
it is also possible to use an interferometric method, and one 
such technique has been developed. This had the advant- 
the entire sheet 


age that thickness variations over the 
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area could be obtained from a single photographic record. 
Comparisons between the volumetric and interferometric 
techniques have shown good agreement, if, in the former 
techniques the coeflicient of velocity is assumed to be 


equal to the coefficient of discharge. 


Application of light interference to the 
sheet thickness 


measurement of 


When a monochromatic light beam is allowed to fall 
on a thin film, reflection takes place from the front and 
back surfaces, (Fig. 8). These reflections interfere with each 
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Warp 
other because of the difference in their light paths, and it 
can be shown that the optical path difference D is given by 


2S COS (11) 


where mm refractive index, 
and x angle of refraction of the light beam in 


the film. 


If D is an odd number of half wavelengths, the two 
reflected beams reinforce each other and a bright fringe 
occurs, Thus from equation (11), the change in thickness 
AS between two adjacent bright fringes is given by : 


A 


AS (12) 
COs 
where A is the wavelength. 
The variation in sheet thickness from a fan spray 


nozzle could thus be determined by measuring the thick- 
ness of the s’. cet in the orifice, i.e. the width of the orifice, 


and the distance of each fringe from it. However equation 


(5) shows tha tnear the orifice, dS/dx is very large, so that 
adjacent fringes will be very close to each other and 
difficult to resolve. This difficulty can be circumvented in 
the following manner. 

In Fig. 9 let the thickness at the rth fringe, (the first 
clearly resolvable fringe) be S, at a distance x, from the 


orifice and let the subseript n denote any subsequent 


fringe. 
Thus from equation (5) 


and 


and therefore S S 


r 
From equation (12) 


2 COS 


Therefore 


n r 


\ linear relationship should thus exist between (7 r) 
and 1 ry, If the slope of the straight line is given by m, 
then 


A 


A 


lm 


Now 


sin 2 sind » where ¢ is the angle of incidence. 


Thus the expression may be rewritten as 
\ 


(14) 
2m pil sin? ¢ p?)! 


Experimental technique for measuring sheet thickness 


The apparatus employed to record the fringe patterns is 
shown in Fig. 10. The light source, M, is a 250 W compact- 
source mercury vapour lamp, whose light is condensed 
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by lens L, and transmitted through a mercury yellow filter 
F and diffusing screen D, onto the liquid sheet. The 
reflected light is photographed by the camera P. Liquid 
is fed to the nozzle N (seen from above) by compressed air. 
Its pressure is measured by a Bourdon gauge and its 
temperature by a mercury thermometer. The refractive 
indices of the liquids have been measured by a refracto- 
meter. The light consists of the two yellow lines of wave- 
lengths 5770 and 5790 A, and the mean value 5783 A has 
been used for calculations. 

Equation (14) shows that the value of K depends on the 
value of sin? ¢. The error in determining 4 can be reduced 
to a minimum by maintaining it at a small value. ¢ can 
be made zero by employing transmitted fringes but 
because these are produced by a double reflection, they 
have a very low contrast and cannot be satisfactorily 
photographed. Reflected light is therefore used and ¢ is 
kept to a minimum by placing the camera and light 
source as close together as possible. Under these conditions 
the value of dé is 11. This is considered satisfactory since 
an error of % only causes an error of O-6 per cent im 
value of K. The camera is aligned at the correct angle by 
a simple parallax method. 

Because the plane of the sheet is disturbed by waves, 
relatively long exposures are required to ensure that light 
is reflected into the camera from every part of the sheet 
as the latter moves past the lens. On the other hand the 
exposure has to be kept to a minimum in order to reduce 
external vibrations. The best compromise is found by 
illuminating the sheet for 1 50 sec and employing Ilford 


Rapid Process Pan emulsion. 


EXPERIMENTAL RESULTS 


The thickness parameter A has been measured 
for one nozzle (Fig. le) over a wide range of 
liquid surface tensions, viscosities, and injection 
pressures. Two levels of approximately constant 
surface tension, 30 and 73 dyn/cm were studied 
by using kerosene—paraflin and water—glycerol 
mixtures respectively. Viscosities were varied 
from 1 to 50cP and injection pressures up to 
50 p.s.i.g. A few determinations were also made 
with ethyl alcohol, benzaldehyde ethyl 
acetate. 

Typical interferograms are presented in Figs. 


13, 22 and 23. Fig. 11 gives a typical plot of 


(n —r) vs 1a, and shows an excellent linear 
correlation as predicted by equation (13). 
Irregularities in the nozzle orifice caused the 
value of A to vary with different streamlines 
(see Plates 23, 24) so that a mean value has been 


calculated for the core of the sheet over a total 


angle of 45° round the nozzle axis. 
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This work has confirmed previous findings [7] 
that the coefficient of discharge of a fan spray 


nozzle is independent of surface tension or 


viscosity. In the present case it has a value of 


Effect of injection pressure on the thickness 
parameter K 


Fig. 12 shows K as a function of injection 
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Fic. 12. The effect of pressure on the thickness 
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pressure for values of viscosity below 2 cP and 
that it asymptotically to 
constant value of approximately 14 w4 
of pressure. Thus, 


nozzle produces a sheet which is less than 


indicates tends a 
with increase liquid this 
thick at a distance of 38cm from the orifice. 
Certain differences can be 
the of high 
(73 dy n em) and kerosene and alcohol, the liqquids 
of 
The curve for water shows that as the injection 
A to 
However, kerosene, 
10 


observe ad between 


liquid surtace tension, 


water, 


low surface tension (25 and 23 dyn cm). 


steadily increases an 


kK 


tom? at 


pressure is raised, 


asymptotic value. with 
rises rapidly to a maximum of 14-7 
5 p.s.i.g. after which a steady decline takes place 
th The 


alcohol appears to be similar to that of kerosene 
although the 


towards asymptote behaviour 


lower myection pressures were not 


examined. 
Liquids of higher viscosity were only examined 


between an injection pressure range of 15-50 


p.s.ig. It was found that the value of A tended 


to diminish slowly with increase of pressure. 


inal ck 
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are ad out. 


Fig. 13 shows a typical photograp! 
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low viscosities (Figs. 22 
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14 shows the 


and at injection pressures of 15, 30 and 50 p.S.1.g. 
respectively. 

Kach curve is plotted for data pertaining to two 
levels of surface tension, and the uniform scatter 
of the 


clearly 


around curves 
that 


mental error there 


their respective 
the of 
is no effect of surtace 
Surtace 


points 


shows within limits expert 
tension, 
This result is to be expected, tension 
should not affect the flow behind the orifice which 
can be considered to be the main factor inthuencing 
the sheet thickness. 

However, the value of A is markedly affected 
14 shows that it 
The rate of 
At 15 
10 the 


\t 50 psig. K 


by viscosity, and Fig. increases 


with increase of viscosity. merease 
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Effect of liquid properties and injection pressure on 
the trajectory of the edge of a fan spray sheet 


Equation (10) can be re-arranged to give 


sin (p whose validity 
can be tested by plotting the parameters ry P.A. 
@ Since the 


origin of the sheet can be considered to be very 


and sin (8 #) since C, 18 constant. 


close to the orifice, 2 has been measured from the 
oritice, 

Fig. 15 illustrates the result for liquids of low 
18cP) and covering a wide 


viscosity (Os to 
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Fic. 13. Interferogram of sheet of viscous liquid (glvcerol-water 50 cP at 15 p.s.i.g.). 


(Exposure time 1, 50 sec). 


a. 
EF 
= 
>! £ 
O-Scm 
. 


4 
4 
VOL 
12 
19¢ 
. 
h 
. 
= 


Theoretical line 


+f) 


The trajectory of the sheet edge at low viscosities 
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N. Domprowski, D. 
range of experimental conditions (injection pres- 
sure 5-480 p.s.i.g., surface tension 25-475 dyn em, 
and density 0-8 to 13-6 g cm*. It shows that the 
experimental points lie within a band whose 
slope is somewhat greater than that of the theore- 
tical line. i.e. the theory predicts smaller sheets 
than those observed in practice. This is to be 
expected since the derivation has not taken into 
account the fact that although most of the liquid 
amassing in the contracting rims is continuously 
being lost as drops are being ejected (cf. Fig. 2), 
the momentum remaining in them will restrain 
the action of the surface tension force. However, 
bearing in mind the complex nature of the prob- 
lem, it can be seen that the basic assumptions 
underlying the theory are justified. 

Fig. 16 illustrates the effect of a viscosity range 
of 1-50cP at a constant injection pressure of 
15 p.s.i.g. and mean surface tension of 70 dyn em, 
It shows the curves to deviate from the “inviscid” 
theoretical line by an extent depending on the 
liquid viscosity. This effect probably results from 
the reduction of the contraction velocity VO by 
viscous forces which were not taken into account 
in deriving the relations. 

Fig. 17 illustrates the effect of liquids at one 
level of high viscosity and two levels of surface 


tension. The figure shows that liquids of higher 
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surface tension deviate further from the theoretical 


line than liquids of lower surface tension. This 
is presumably a consequence of the greater 
viscous retarding forces operating at the higher 
contraction velocities for liquids of hizher surface 
tension, 

Fig. 18 illustrates the effect of viscosity at 
different levels of injection pressure and shows 
that the theoretical curvature is approached as 
the pressure is increased. This results from the 
fact that at high injection pressures the high 
radial velocities limit the extent to which the 
sheet can contract. Deviations from the th-ore- 
tical contraction velocity will therefore only have 
a slight effect. 


Discussion or RESULTS 


The effect of injection pressure, surface tension 
and viscosity on the thickness of the sheet can be 
attributed to two distinct causes, namely the 
conditions of flow inside and outside the orifice. 

The ftlow imside cireular orifice can be 
characterized by the Reynolds number based on 
the orifice diameter [15], and therefore under 
conditions where the orifice tlow controls the 
flow in the sheet, the thickness, as expressed by 
AK, should be capable of correlation with the 
Reynolds number, 
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Fic. 17. The effect of surface tension on the trajectory of the sheet edge at high viscosity. 
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Fic. 18. The effect of ejection pressure on the trajectory of the sheet edge at high \ iscosity. 


(Viscosity 50¢P, surface tension 70 dynes cm) 
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The present series of experiments were carried (WV). can be grouped together in the form 


out with only one nozzle size so that the effect Vp 7 RO n since the coetlicient of discharge 
of the linear dimension in the Reynolds number is constant the parameter R can be written as 
cannot be verified. However. the remaining (pp)? », where P is the injection pressure, 

factors of Viscosity, (»), density, (p) and velocity This group has been determined for all the 
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experimental conditions and ts plotted on Fig. 19. has been confirmed from values of sheet thickness 
The results show the points to lie satisfactorily reported earlier [4] for water and a series of 
around a smooth curve except the values for water nozzles of the same design as that used in the 
and kerosene at low injection pressures (cf. Fig. 12). present work. These data, converted to values 
As the magnitude of R increases A tends asym- of A. are plotted in Fig. 20. The angle, @,, 
ptotically to a constant value of about 14 10-4 caleulated from the slope is 50°. The value of @, 
em’. calculated by equation (6) for the nozzle used in 
K is also related to the angle 0, by equation (6), the present work is plotted in Fig. 21 as a functtion 
20, being the total sheet angle if the edges were — of (Pp, y)"* and it can be seen that G7 asynipto 
unconstrained by surface tension and the orifice tically approaches a maximum value of 75°. The 
exit passage. The equation predicts a linear exit half angle of the clearance passages (see Fig. 
relation between the orifice area and A. This le) is 60°, so that for the larger values of 0+. the 
edges of the sheet must be constrained in order to 
leave at the former angle. This IS contirmed by 
measurements of the sheet angle at the nozzk 
exit. The inconsistency in values of the two 
angles, for nozzles of nominally identical design, 
im the present and previous work must be due to 
minor differences in the flow passages. This could 
not be easily ascertained because of the small 

size of the specimens, 
The observed variation of A and @, with flow 
conditions can be attributed to the movement of 
4 . the centre of pressure behind the orifice. At a 
particular value of R. the centre of pressure ts 


fie. 20. Relation between orifice arca and parameter K situated in such a position that the licquaicl flows 
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directly along the exit walls (point {. Fig. 1). 
Under these conditions the value of the angle Oy 
is identical with the angle @ (Fig. 5) at which the 
sheet leaves the orifice. At lower values of R, the 
entre of pressure moves upstream since viscous 
Heag reduces the effective angle of convergence. 
Thas results in a thicker sheet (i.e. a larger value 
of At) diverging at a relatively small angle. An 
incre ased value of R causes the liquid streams to 
impitage at a larger angle so that the sheet tends 
to diverge to a greater extent as the centre of 
pressure moves towards the orifice. However, 
although the value of the angle @> is increased, 
the sheet edges are constrained by the exit walls 
so that the actual angle of the emerging sheet 
can Increase further. 

‘The results for kerosene and water at low 
injection pressures (Fig. 12) indicate an unusual 
flow pattern which cannot be explained by the 
conditions in the ortlice and must be due to an 
external cause. possible explanation lies im 
the fact that at least low pressures the spray 
sheet is complete ly bounded by thick rims. Under 
these conditions, surface tension § forces are 
comparable to the momentum forces of the sheet 
and it is likely that the normal flow pattern in 
the sheet is affected. 

This hypothe SIS IS partly confirmed by examina 
tion of the interferograms of kerosene and water 
in Figs. 22 and 23 respectively. At low pressure, 
kerosene (surface tension 25 dyn/cm) exhibits 
fringe contours which are concave with respect 
to the orifice. As the pressure ts increased to 
5 p.S.Lg., the contours flatten out and then 
reverse to become convex. Since each fringe 
represents a line of equal sheet thickness, the 
photographs show that below 5 p.s.i.g., the 
central sector of the sheet is thinner than the 
surrounding parts and that the flow only 
approaches the ideal case at higher pressures. 
On the other hand, the fringe contours for water 
(surface tension 73 dyn cm), corresponding to the 
asymptotic increase of A shown in Fig. 12, change 
their character only gradually, and remain 
concave at the highest pressure. 

The occurrence of this flow pattern at low 
injection pressures probably accounts for the 


observation that at verv low values the thick 


rims occasionally break away from each other to 
form two liquuid jets. This is also a possible 
explanation for the fact that with liquids of low 


density and high surface tension such as sodium, 


the sheet is disrupted almost as soon as it is 
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Abstract 


of phosphorescent slurry particles into the reaction medium. 


Certain radiation promoted chemical reactions may be accelerated by the introduction 
These particles adsorb radiation 
energy in the regions of high radiation inte nsity and release this energy throughout the reactor. 
A solution for this distributed energy is obtained in terms of the properties of a stirred reactor 
radiation system and the results applied to a design problem. 


Résumé 


lérées par lintroduction de particules d'une boue phosphorescente dans le milieu réactionnel 


Les réactions chimiques activées par une radiation quelconque, peuvent étre accé- 


Ces particules adsorbent l énergie de radiation dans les régions d’intensité de radiation élevée et 
libérent cette énergie dans le réacteur. Une solution de cette re partition de lenervie est obtenue 


en fonction des propriétés du systéme de radiation du réacteur A agitation ; les résultats sont 


appliqués au caleul du réacteur 
Zusammenfassung —Kinive durch 


Strahluny Reaktionen 


durch Kinftihrung phosphoroszierender Teilchen in das Reaktionsme dium beschleunigt 
hese 


hervorgerufene chemisehe 


kOnnen 
werden. 


Teilchen absorbicren Strahlungsenergie im Bereich hoher Intensitit und strahlen diese 


Energie innerhalb des Reaktors aus. Die Energieverteilung wurde berechnet fir den Fall eines 


yverthrten Strahlunygsreaktors und die Ergebnisse fiir ein Entwurfsproblem angewandat. 


INTRODUCTION This paper presents a method for the caleula- 


} 


tion of the average distributed energy 
of a series of continuous irradiated stirred 
tank reactors and extends the method to the case 


has been increased interest In eac 


of 


radiation. 


Recent.y there 
the 


energy 


in 
chemical unit 
In 
large portion of the reaction medium is shaded 
the 


become activated in the regions of high intensity 


promotion reactions by high 


industrial size reactors a 


of a tubular reactor with axial dispersion. 


from radiation source. Slurry particles 


THEORY 


radiation and then lose this activation ene rey to 


many types of reactions the effectiveness of the 


In irradiated continuous stirred tank reactors, 


the chemical reaction in the shadow regions. slurry particles become active in regions of high 


radiation and then lose that 


intensity 


energy 


radiation depends on this distributed en rev. 

An interesting example is the proposed scheme 
for the purification of space vehicle or submarine 
The 


carbon dioxide laden atmosphere is fed into a 


atmospheres by a photosynthesis reactor. 


series of chemical reactors in which light radiation 
promotes a photosynthesis reaction producing 
oxygen. Phosphorescent particles are introduced 
into the reactor to increase the efliciency of the 
light. 


light source and then re-emit their energy in the 


These particles become active near the 


shadow region. 


*Present Address 
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by decay to the chemical reaction. The complexity 
of this activation process suggests a probabilistic 
the the 
distributed energy. 
Let be the 
particle at time € if it was activated at time zero. 
This the at the 
particle loses energy and is obtained by measure- 


approach to solution for average 


g (€) energy associated with a 


function describes rate which 

ments on single isolated active particles, 
The radiation is attenuated by the reaction 

medium: those particles nearer the radiation 


source have a higher probability of being activated. 
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Let P denote the probability per unit time that a 


particle will become activated, This parameter 


depends on the intensity of the source, the 


mixing properties of the reactor and the opaque 
medium. In the case of 


ness of the reaction 


perfect mixing it can be estimated without 
difhieulty. 
Define p(r) da as the probability of a party le 


bemg cuM of 


at position rtoa dy. For the 
pertect mixing and axially symmetric geometry 
p («) da dv H where H is the depth of the 
reactor. 

f(a) is the that a 

This 
proportional to the 
the 


SCOPIC 


probability per unit tonne 


part le at will be activated 


position 


probability function — ts 


radiation intensity at position v7, pro 


portionality constant being the 


absorption cross section, 8. of the particles, 


Bi ir) 


party le ws al position 


rtoa 


probability that 


dx and is activated in unit time is then 


{ (1) p(a) da 


integral of this function over the entire 


The 
reactor is defined as P. 


H 


BI (a) de. 


for perfect mixing and if the radiation ts absorbed 


according to Beers law, 


I (x) = I, exp (— 8), 


P =I, 8/s{[1 


exp ( 


where J, is the radiation intensity at the source 
and s is the linear absorption coefficient of the 
reaction medium. 

The age of a particle ts defined as the time 
since its last activation. This distribution of ages 
can best be described by a distribution function 
the 
The 


average activity per unit mass of particles is 


f(é). is defined such that EV dé is 
fraction of particles of age to 


then the weighted average 


distribution of ages in the k"" reactor of, 


The 


the series of determine 


tors mitist now be 
in terms of the properties of the reactor-radiatAon 


This can be 


a material bal 
The 
the particl s activated in an merement of tine 
dr is Pdr and the fraction that escape activation 
Pdr. 

The total mass of particles of age & to ¢ dé is 
f, da where is the total mass of ticles 


svstem. dome by 


over the age to fractigon of 


yroup < 


is 


reactor and the k refers to the 
“tate 


inn thee 


reactor the steady 


series. In the 
the total mass of particles in the age yroup that 
appears in an increment of time dr must equal 
the mass that disappears. 

Particles in the age group € to ¢ d& way 
appear by : 


1. Natural 


amount (1 


avemg of particles mn an 
Pdr) drjd 
into the reactor from the 


dr f, da where w 1s 


reactor mm an amount 


the mass flow rate of particles, 
The particles disappear by ; 
1. Activation in an amount Pd Wy, (£) dé, 
2. Natural Pdr) 
Wf, (€) dé. 
3. Flow 
wd + f, (€) dé. 


ageing im an amount (1 


from the reactor m an amount 


mathematical transliteration of — the 


The 
material balance ts 


dr) dr 
(&) 


(€) 


(1 Pdr) Wf, 
wadr 


which becomes upon rearranging 


f, (&) Fats dr) 


dr 


Pf, dr) 


Kquation (2) becomes as dr > 0 


(£) 


dé (8) 


(P + (€) 


| 
G = (1) 
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where = w W. 


Particles of age 0 to 0 d§ appear by - 


1. Activation of any particle in an amount 
W. 
2. Flow into the reactor m an amount 


wd +f (0) d€. 


1 


The \ disappear by ; 


1. Natural ageing in an amount 
(1 Pd +r) W f, (0) dé. 
2. Flow out of the reactor in an amount 


wd rf, (0) dé. 
The material balance on this ave vroup becomes 


as dr -Qand di 
(0) P 1) 


Equation (4) is the boundary condition necess- 
ary for the solution of equation (3). Equation (3) 
is a Standard linear ordinary differential equation 


hose solution Is 


f, (© = Pexp(—[P + Q)& 


exp | (P 


i(é n)dn (5) 


Equation (5) shows how the age distribution fune- 
tion changes from reactor to reactor. This allows 
the calculation of the age distribution in the 
reactor in terms of the properties of preceding 
reactors and the age distribution in the feed to 
the first reactor. The mean activity in the k™ 
reactor is obtained by use of equation (1). 

These results can be extended to the case of an 
irradiated tubular reactor with axial dispersion. 
It has been shown [1] that a tubular reactor with 
axial dispersion can be approximated by a mixing 
cell (stirred tank) model. The length of the 
mixing cells must be chosen so that uAv/D 2 
where Av is the size of the mixing cell, u is the 
average velocity in the reactor and D is the 
axial dispersion coetlicient. The average activity 
at position Lin the reactor ts approximately 
that given by equations (1) and (5) where & ts 


taken as L Av, 


Energy distributions in irradiated slurry chemical reactors 


An EXAMPLE 

The application of these theoretical develop- 
ments can best be seen in a design problem. 
Consider the promotion of the biological reaction 
{ +B in a series of illuminated continuous 
stirred tank reactors. The primary illumination 
activates phosphorescent particles which re-emit 
secondary radiation. This secondary radiation is 
of the proper wave length to promote the chemical 
reaction. How does the fractional conversion of 
this reaction system depend on the operating 
parameters ? 

The active particles lose their energy by eXx- 
ponential decay : that is, g(&) G exp | ~&) 
where G is the energy absorbed by the particles 
at the time of activation. 

kxperimental considerations show that the 
rate of reaction of specie A is of the form 


n 


Rate = C, Rp 


where C , is the concentration of A, p is the mass 
of particles per unit volume of reaction medium 
and both R and w are constants determined 
experimentally. 

If the feed to the first reactor is inactive, 


equation (5) gives the age distribution in the i 


reactor as 


k 
VY 
f, (€) = Pexp(— €[P > 


Using equation (1) the average activity per unit 


mass of particles in the k'® reactor is 


\ material balance on the chemical specie A 


about the &“" reactor gives 


1 
R[p(G, 


tk | 
( i i 
where 4 is the nominal holding time of the 
reaction medium. 

Equation (7) gives the concentration in the i 
reactor in terms of the feed concentration C, 


and the properties of the preceding reactors. 
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Rupp 


PG, (W1,8/ VHs) (1 — exp — 


G Q 8/Hs) (1 exp sil) 

The fractional conversion 1—-C/C, 
is the fraction of the reactant A that has been iF) ™ 
converted at the k™ reactor. Using the relation (1,8 Hs\(1— exp — sil) 4 | 
for P as well as equations (6) and (7) the fract- ; 
ional conversion at the k" stage is given by 


equations (8) and (9). 


R [p (Gs/G)}"| 8) 


where 


mun 


Pia. 3. The effect of the flow rate of particles on the 


fractional conversion 


The effects of the various operating parameters 
can best be shown graphically. Fig. 1 shows the 


Fic. 1. The effect of the mass of particles per unit velume — Cffeet on the fractional conversion of the mass of 


retained in the reactor on the fractional conversion ba 
particles retained in each reactor. Fig. 2 shows 


the distribution of ages in each reactor in the 


series. The effect of the mass tlow of inactive 
particles is shown in Fig. 3. 

The effects of many other parameters can be 
studied (e.g., source intensity J5, depth of the 
reactor H, etc.). These studies give insight into 
the performance of such a reactor system. 

The physical constants used in the example 
are : 

100 
O-285 
OS 
I, B/s [1 exp | H = 0-2 min~! 


3-5 min 


The distribution of particle ages, a == 1-0 min"! 
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CONCLUSION model provides the insight for the experimental 
A mathematical model has been developed work that must follow. 
which isolates the effects of the important 


icknowledgements—1 wish to thank STerHen 
operating parameters in a series of irradiated 


and Ruruerrorp Artis for their suggestions on the 
Apntinuous stirred tank slurry reactors. This mathematical presentation. 


REFERENCE 


[1] Coste J., Rupp D. F. and Amunpson N, Taylor Diffusion in Tubular Reactors. In press. 


| 
iy 
VOL. 
960 
55 
4 


Chemical Engineering Science, 1960, Vol. 12, pp. 56 


Pergamon Press Lid., London. Printed in Great Britain 


On Denbigh’s optimum temperature sequence 
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Abstract The example of an optimum temperature sequence in stirred tank reactors given by 


Denpicn is here discussed from the dvnamix programming pomt of view The notion of dvnami 


programming was developed t leal with just such optimum problems, in which a decision 


in this case a choice of te mperature, has to be made at each of N stages that will Optimize scone 


output of the whole Process \ general presentation of the principle is given but app licat ie nm is 


made only to Dennian’s example. The solution is extended to the case where physical restrictions 


ire tnposed on temperatures and holding tines 


Résumé loexemple donne dune sequence de température optima dans des 


reacteurs a curve comportant une agitation est 


dliseute ici du point de vue de létablissement 


de la programmation dyvnaniqu Cette notion est developpee dans le but de traiter exactement 


ces problemes ou il faut decider du choix de la te mperature a chacun des N étages cette tempera 


ture doit fournir le meilleur rendement pour Popeération global Une presentation générale du 


principe est donnee mais lapplication n'est faite 


que sur Texemple de Denier. La solution 


du probleme est étendue au cas of des restrictions | hvsiques sont imposées aux températures et 


aux temps de séjour. 


Zusammenfassung Das Beispicl einer optimaten Temperaturfolve in’ Rihrwerksbehiltern, 


die von entwickelt wurde, wird vom Standpunkt der dvnamischen Programmicrung 


diskutiert. Die Kone zption der dynamischen Programmierung wurde zur LOsung von 


rungs aufgaben entwickelt, bei denen cine Entscheidung (in diesem kalle die Wahl der T mperatur) 


auf jeder der N Stufen zu treffen ist, um die Ausbeute des ganzen Prozesses zu « ptimieren. Eine 


allgemeine Darlegung des Vrinzips wird gegeben, die Anwendung erstrecket sich jedoch mur 


auf das Deneiu'sche Beispicel. Die Losung wird fiir den Fall erweitert. wo phvsiaklische Kin 


schriinkungen fiir Temperaturen und Verweilzeiten vetroffen werden 


INTRODUCTION 


Ir HAs been rece ntly shown by Denpicn that 
remarkable improvement in yicld may be mad 
by corres tly choosing the temperatures at which 
a sequence of stirred tank reactors should operate. 


The system of reactions he considered was 


(4) 


in which is an intermediate. VY the desired 


product and P and @ waste products. It is 


assumed that 4 is a costly raw material so that 


it is of importance to minimize the production of 


P and ( or maximize the vield of Y as a fraction 


of A. Uf the reaction (1) has a lower t mperature 
coetlicient than (2) it will be desirable to operate 
at a low temperature before any appreciable 
amount of \V has formed, but if reaction (3) has a 
hicher temperature coetlicient than (4) a high 
tompe rature will be desirable in the later stages 
of the reaction. This suggests that two reactors, 
the first held at a lower temperature than the 
second, should give an mnprovement of the vi ld. 
low vreat an improvement this may be is shown 
by Dennicn’s example, in which a single reactor 
at 360 °K gives a 25 per cent vield, whilst two 
reactors the first at 280 °K and the second at as 
high a temperature as possible give nearly 
58 per cent conversion of 4A to Y. It was assumed 
that the second reactor was so large that the 


reactions, all of them first order and irreversible, 
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go virtually to completion and the holding time 
of the first reactor was taken to be the time at 
which the concentration of VY in a batch reaction 
would be greatest. Denpicnu’s method was one of 
clirect calculation for various values of T, and 
the temperatures in the two reactors. He re- 
cognized that the holding time in the first reactor 
should also be varied but was able to show that 
the yield was not very sensitive to this. Even 
with two reactors, the method of direct calculation 
might, with less perspicacious handling or 
constants not so nicely chosen, founder in a morass 
of interminable approximation; if more than 
two reactors were considered, more than ordinary 
insight into the behaviour of the system would 
bn required lo organize such a cal ulation. Some 
interest therefore attaches to the probe moot re- 
ducing thes optimum calculations to a simple 
algorithm such as could be performed by a com 
puter, human or mechanical, with no insight into 
the phy sical probabilities. 

An organizing principle ideally suited to such 
problems is provided by the notion of dynamic 
This has 


and his associates | 2] and applied to a 


programming, been developed by 
wide variety of multi-stage decision problems. 
This is just the kind of problem that confronts 
us; each reactor is a stage and at each stage a 
decision as to size and temperature has to be 
made, which will optimize the final yield. In 
the next paragraph we will discuss the relevant 
part of dynamic programming in gene ral terms, 
but return to Dennicu’s example to show how 
It will be obvious that the 


provides method of solving almost all 


it works in practice. 
theory 
optimum problems of this type. We shall not 
discuss the general problem here, for much of 
the interest and most of the instruction in any 
theory lies in specific cases, and with so elegant 
an example to hand it would be stupid to over- 
look it, 


problems in reactor design will be pive;n ¢ lsewhere, 


An extended treatment of optimum 


Tue Norion or Dynamic PROGRAMMING 


Consider a process consisting in a finite number 
of discrete stages N. The state of the process 
described 


(Par Paz Pum) 


th 
stream leaving the xv‘ stage can be 


by a set of m quantities p,, 


and the feed to the first stage is in the state py. 


" stage be 


Then the effect of each stage is to transform the 


Let the operating conditions at the vn‘ 
given by a set of r quantities q,, 


state of the incoming stream, p,, into the state 


P,» this transformation depending on the operat- 
ing variables g,. Accordingly we may write 


Pn T (p,, Vn) (Py, 1) (1) 
In the present context we may think of p, as 


concentrations of components, g, as holding 
time and temperature and the transformation 
reaction, The 
Po + + + Pn 


po by a series of transformations 7 


as the sequence of states Pp 
is generated from the initial state 
f,...f 
which «ce pe nd on the choice of ope rating variables 


Un 
bles ts called a policy, and if the g, satisfy certain 


Any choice of the operating varia 
imposed restrictions (e.g. temperatures must be 
with certain limits) the policy is called admissible. 

We now detine a profit function | P(py) — P(p,)! 
which is a measure of our satisfaction over the 
change of state resulting from the process ; in the 
present case this might be the yield of the valuable 
product. The admissible operating policy which 
MAXIMIZeS P (py) for given initial state is said 
to be the optimum N stage policy with respect 
to the 


bonum we seek. When this policy has been found 


initial state py, and this is the summum 


the resulting maximum of P is a function only of 
the initial state p, and we may write 
Max } P (py) 


fy P (2) 


the maximization being over all admissible 


polien If we were to seek this maxir m 
directly we should have to vary all the rN ponicy 
variables simultaneously and this would be a 
task. 
the N stage process as a first stage followed by 


an 


the choice of 


considerable However, if we consider 


1) stage process, It is clear that whatever 
operating variables in the first 


stage, the overall policy cannot possibly be 


optimum unless the remaining (N 1) stages use 
the optimum (N 1) stage policy with respect 
to the state p, resulting from the first stage. 
Moreover by varying 4, and so py we shall tind 


the V1 which together with the optimum (N 1) 


: 
| 
: 
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R, 
stage policy constitutes an optimum N° stage 
policy 

Since P (py) 
In an equation 


P (po) ={ P(p,) — P(p»)} 


P (p,) {| we can express this principle 


Max (p,) 


Iw (Po) P + (3) 


where 


Py (Po) (4) 
the 


admissible policies q\- 


the 
Equation (3) may be 


and maximization is now only over 


transliterated by saying that 


‘The maximum 


profit from N in 


Maximum 


reactors with 


\ feed Vo 

where p, is the state of the product from the first 
reactor with feed p, and operating conditions q, 
and the maximum is sought by varying g,. The 
principle of optimality has reduced the number 
variables that have to be simultaneously varied in 
rN to Since the 
amount of work depends exponentially on the 
of 
this is a most important saving. The sequence of 
the 


seeking a maximum from 


number simultaneously varying quantities 


equations (3) must be started by solving 


single optimum 


(Me) = Max { P(p,) — P (py)! (5) 


In seeking an optimum policy with respect to 
had 
imbed the problem in a rather larger one, for at 
rate until the (NV 
determine f, (p) for a whole set of 
of Py 
domain of The 
problem is still far from trivial but it has been 


the particular feed state Po we have to 


any 1) stage we have to 
states p, 
general variation will make 


Since 


vary over certain p-space. 
reduced to a straightforward algorithm and one 
well-suited to repeated application. In addition 
the functions and policies generated on the way 
are not without interest or use, for it would be 
a short-sighted analysis that considered only 
one initial state py. Furthermore restrictions such 
as upper and lower bounds on the variables q,,, 
which are diflicult to handle analytically, actually 


make the problem easier from the point of view of 


‘The profit made 
the 


actor with feed 


ARIs 


mechanical computation, for they confine the 


area of search. 


3. Tur Equations 


We will denote by a, x and y the concentre.- 
tions of A, X and Y and let ¢ be the holding 
time of a reactor. If suflix n denotes the values 
at the n™ stage a mass balance gives the three 


equations 


t, (Ay, kes, ) ay (6) 


t, ky, a (ks, Kan) (7) 


n 


The maximum profit 

that 

by (N 
with feed p, 


first re- can be made 


1) reactors 


(8) 


k,(T.,) is the rate constant of the i 
th 


Here k,,, 
reaction at the temperature of the n»" reactor. 
To reduce the number of rate constaats we follow 
DeNbIGH in using the ratios of the reaction rates. 
Let 

r, /ky tk 


In 


(9) 


Since r,, is a function of 7, only and is monotonic 


in 
we can let it play the role of temperature in what 
follows ; s 
is known. 
Specifically, if 
T 


, Will then define the holding time when 


k* exp — E,/RT, 
E,)/ R ln (k,* (10) 
and 


(11) 


“a i'n 


The variable r is the ratio r, in Densicu’s notation 


and the other ratios rz = k,/k, and ry ky/k, can 
be expressed in terms of r in the form r; 
(E, E,)/( FE, E,) 
(E, — E,)k,* +( 2, — E,)k,* 


If we were to do a complete 


P; 

where, 

and (FE, — E,)inp, 
(E, — E,)k,*. 


study of all reaction systems of this type we 


for example 


should have essentially four parameters pg, 75. 


p, to vary. This might be carried through the 


analysis that follows but for simplicity we will 
immediately adopt Dennicn’s values 


rs = 107? (12) 


| 
= 4, + 
\ \ q 
| —— VOL 
12 
Po / ] 
Yn Yn Ks, rn 
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Equations (6)(8) now become The variable a,, is related to the fraction of A 
which is used in the n™ stage and it simplifies the 
algebra slightly to use it in place of s,. 4, 1S 
thus playing the role of the holding time ¢, and 
this is perfectly satisfactory since the Jacobian 
d(T,,. t,)/2 %,) nowhere vanishes. 


As a final simplification we notice that these 
Before formulating the dynamic programme we 


equations are homogeneous of the first degree so 
that we can deal in the ratios of the concentra- will consider the optimum conditions in a single 
reactor. We wish to maximize m,, i.e. the yield of 


tions. 
Y as a fraction of A. Firstly it is easy to show by 


Let 
differentiating that », is a monotonic increasing 


function of «, so that , should be made as large 

(13) as possible ; in fact, the assumption that the 

(y, Ya-y)/8n-1 Dp reaction goes to completion is equivalent to 
ty > 

VOL. then the preceding equations may be solved to Secondly by differentiation with respect to r, 

Le Rive and rearrangement we find that for any given 


ws (14) 7, IS Maximum when satisfies the equation 


200 a, (1 
(18) 
(% 100 r,*) (1 + 7) 


where 


L 


, For constant &, the solution of this equation Is a 
" l+r, curve in the «, r plane lying between the extremes 


for &, = 0 and o, namely 


200 72 (1 ry) 


and 
x,,)(1 100 100 r,* 


Optimum conditions for a single reactor, 


En L,, 1 M, (15) 
1 
¢ 
4 So j 
"In 14, (16) l + 
a 
M,—(14 (1 (17) 
1 
a 0-5 4 é | 
g | 
| 
| 
2 


R. 


Unfortunately these lie in a rather narrow band 


even m the low x. log 7 plane (sex Fig. 4) and to 


show bn tte r Fig. l Is 
1 the 


have vertical asymptotes at r, (1 


how thev run a plot of 


low (x 100 r?) versus log r. For curves 


for &, | they have horizontal asymptotes. The 


optimum temperature for any x, and & is given 
by the value of r, at the intersection of the curve 
with appropriate 


ot the 


and the straight line 


ome lines is shown. 
When the reaction goes to comple thon we can 
obtain « xphieit formulae, 


(14); 


for as 


16) and (17) then vive 


whence the conditions 


infinite, its limit being 


shall 
that it is 


cannot be 

Wi 
as in example, 
that the 


necessary 


Actually =r, 


as 7 > 


> however, 


very large 


formula (22) may be used. 


The 
are not diflicult 


simple 
modifications to make (22 


to make. 


exact 


Dyn Pron RAMMING Fort LATION 


In fact r, cannot be infinite as the 


pohiey demands, 
T +x. We shall 
example, that this ratio is so large that r, (1 ’ 
to 1 


is verv small, 


optimum 


k,* as 


for it lie S 


assume, as Deneicnu’s 


} 


is very close for high temperatures and 


r, (1 r; Then we can retain as 
formula 


to 


a very good approximation the 
the 


simple 
modifications 
to make. 


m (22); necessary 


make not difficult 


Let 


this exact 


In (4. 9) Max (YW Yo) 


where the maximization is overall admissible 


policies ff.) N. Since a 


» is 
the feed composition to the second stage and this 


is the first of the (.N 1) stage process as described 


ARIs 


in Section 2, the principle of optimality given there 


asserts that 


fy (Gq. %) Max fy, — ty + (24) 


the 
lirst stage (2), 
(24) that the 


maximum 
reactors is the maximum of 


where maximization is now over admissibt 
In words equation 


from N 


‘what can be made 


states vield 
in the first reactor plus the maximum vield from 
the feed produced by the first reactor.” 

We can use the homogen¢ ity of the equations 


to simplify this equation by writing 


lg By (So) Vo). 


then dividing through (24) by a, 


ives 


Yo a 
a) 


Max 


1\>1 


Max 


i ~ 
l+a 


\ further modification is quite valuable before 


setting down the final equations. Since a, is 


continually decreasing the ratio & may becom 


quite large, and this will cause dilliculty in 


If we 


the reciprocal of & we can use it when ¢ 


interpolating for also introduce 


1 and 


so always keep the arguments of the functions in 


the range (0, 1). 


Let 


and a (26) 


a 


(24) have four 


used according as & and & 


>1 


then in place of equation we 


equations to are 


reater or less than unity. 


= 
| 
4 
- 
— 
; « 
(1 l ry 
VOL. 
ly ly 
" a 
= (V/A) + (22) 
c 
= Se (25) 
= Z, (€) Max ,(é)! 
= <1. & > 1: 
ay } 
(27) 
| So | 
(Co) Max So 7 1 
a 
h. M: Si se | 
} 
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The terms other than g or # are calculated from 
calculated from equations (15) and (16) and all 
maximizations are over variation of a, and r,. 
complet reaction, 


For a single tank and 


x ©, equations (22) give 


(1 4)(1 hy (Cy) (28) 


For N we 


expressions and the 


cannot to find analytical 


equations (27 


must be 


solved by numerical computation. Figs. 2) and 


sum. arize the results under Den Bicu’s assump- 
tion of complete reaction in the last reactor 


a, — O In Fig. 2 the functions g,. 4,, g, and h, 


are shown and the single point g, (0). The values 


of g, (0) 0-25 and g, (0) O5S7T4 contirm the 


hia. 2. vield for N 


Table 1. 


( omparison 


On Denbigh's optimum temperature sequence 


values calculated by Denpicu and the fact that 
g, (0) 0-663 shows that the three reactor system 
is a significant improvement on the two. The 
optimum policies for one, two and three reactors 


are compared in the following table. 


/ 
VJ 


Fic. 38. Optimum policy for N 


The value soty have been convert al tole rature 
E,) and (k,* k,*) from 


values to ¢, it Is 


using the constants 
Denpicn’s paper. ‘To aseribe 
necessary to assume specific values for E, and k,* 
and use the optimum temperatures in equation 
(11). 

Fig. 3 shows how the policy depends on & or Co. 
In the upper section r, is shown, this is the first 
reactor temperature in a sequence of two reactors, 
Part of the curve for N 1 is shown but no 
attempt has been made to include much of it as 
it has a simple formula, (22). The middle section 
shows the policy with regard to s, with two 
reactors; 8, is infinite by the assumption of 
complete reaction. 
The lower section of Fig. 3 gives the outlet 


composition ratio (, from the first tank under 


of optimum policies 


Te m pe ratures 


times 


A 
; 
0-5} 
for 
61 


the optimum policy. This last graph is needed to 
use these results in the three stage policy. For 
supposing (27) to have been solved for N = 3 
and say £ = 0, we have the optimum values of 
r, and a, and hence of s, for the first reactor from 
the calculation but need to know what the policy 
for the other two reactors should be. The calcula- 
tion has given the value of &, or (, at optimum 
conditions (in this case {, = 0-84) and this is the 
inlet to the remaining two reactors which employ 
optimum policy. We therefore treat 0-84 as the 
{, for N 2 and get for the second r, — 0-51 
(point A), 8» 3-4 (point B), Le. the second 
reactor must be treated as the first of an optimum 
two stage reactor. The lower graph gives ¢, 
(O-84) 0-1 (point C) and this is the outlet from 
the first stage of a two stage system, that is, the 
inlet to the third reactor. This we treat as (, the 
inlet to a one stage reactor which must be working 
under optimum conditions, or since it corres- 


ponds to & 10 1 the third reactor must 


have r, = 8, D. 


Ture Orrimum Prosiem Unper 


RESTRICTIONS 


It is of interest to see how the policy is changed 
and the method of getting it works out under 
some obvious physical restrictions. There will 
rtainly bn upper and lower bounds on the 
temperature, say 7’, 7, T*, and an upper 
bound on the holding time, i, {*, since it ts 
not practical to merease the holding time in- 
definitely by lowering the flow rate. Thes« 


restrictions imply. 


(Zo) 
and 


where q E. { E,) 


i 


hence 2, <s* (1 (30) 


These restrictions define an area in the ar plane 


as shown in Fig. 4, since by the assumptions 


q 0. Also shown on this Fig. are the curves 
£, = 0 and & = o from equation (18), between 
which all optimum conditions for the single tank 


Rh. Arts 


Fic. 4. Nature of the restrictions on the policy, 


lie. Indeed since we know that x, should be made 
as large as possible the optimum conditions for 
given & must lie at the intersection of the curve 
corresponding to the given & and the upper part 
of the boundary, 4BC. As drawn here there is a 
value €, of & for which the curve goes through 
the point B. For & < &, the curves will intersect 
the are 4B, meaning that the restriction on holding 
time is the controlling one. For &, > &, the 
curves intersect on BC, where the temperature is 
the controlling restriction. 

To show how this works out in Deneicn’s 
example we have taken q =5 
and s* $20, for which Fig. 4 shows the nature 
of the restrictions : E, = 0-914. To calculate 
g, and h, it is best to work backwards from a 
given point of the are ABC, calculating £, from 
equation (18) and », from equation (16). These 
functions are shown in Fig. 5 and it is at once 
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On Denbigh’s optimum temperature sequence 


Optimum vield with restricted policy. 


evident that the restriction on temperature 


(T* 394) has an appreciable effect on the yield. 
To calculate g, and h, equations (27) were again 
solved numerically, this time with the maximiza- 
tions over all r, and «, satisfying the restrictions 
(29) and (30). Fig. 6 shows the policy for the 
restricted optimum problem and may be com- 
pared directly with Fig. 3. The following table 
summarizes the optimum policies and yields for 
the restricted problem with & — 0, and the values 
of r have again been translated into temperature 
using Denpicn’s constants. Since we have 
assumed q 1 we can translate the values of s 


into the ratio of holding times of reactors, it is 


Optimum policy with restrictions, 


of course not possible to get absolute magnitudes 
until A,* is given. Thus in the three tank 
optimum the ratios ¢,:t,: ty 3-4: 100: 1. 


. 


6. Disc SSLON 


By modifying the profit function we might 
solve a number of similar and perhaps more 
realistic optimum problems. For example if we 


Table 2. Comparison of restricted optimum policies 


Temperatures 


* Hlolding times 
Yield 


| Aig 
| 
| 
1O} 
| 
~~ 
N2 | 
OF 
| 
| 
8 
0 0-5 
A 
5, 
¥ 
. 
or 
. 
| 
1 | 318 340 22:1 
| 
2 276 2-82 2100 45-1 
3 260 28S 394 1-4 1-3 2100 19-5 
4 
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associate a cost, c, with 4 anda value, ov, with the optimum vield from N reactors 


product Yo owe might seek the maximum of ty (49. 


‘ (Yn Yq) C (le ay). whic h Is lite rally In rg) 


reaction rate reaction 


exp | RT) 


functions defined by equation (17) 


the net profit of the reaction. Alternatively it 
would be possible to allow for the cost of the 


reactor system by taking it to be some function of 
. number of stages 


the holding times ¢.. In this case the formulation 


we have given would not be the most suitable and profit function 


it would be better to work directly with the state of stream leaving the na” stage 


holding time and temperature, The particular operating variables at the nO Stage 
form we have chosen is the simplest for hand yas constant 


calculation, for machine calculation it might well Ky hey 


be better to work from equations (tj S) without upper bound ons 
transformation. However, tt is hoped that lower hound on 
sullic rent has been sand to show how \ 
tomperature 


tool the notion of dvnamic programming ts in the 
holding 


study of such optimum problems. 

. I I upper bound on holding time 
icknocledgements am indebted to Mr concentration of X 
Lee of this university for help with the calculations in concentration of Y 


thos puaper cletimed ty equation (V7 


NOTATION 


concentration of 1 


wetivation energy of reaction 
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Shorter Communications 


Vapour liquid equilibria for ethanol—water-—ethylene glycol system 


( Received 10 May, 1959) 


Abstract——Vapour- liquid equilibrium data for the system ethanol water were determined at one 


atmosphere pressure in the presence of ethylene glycol using an Othmer still in order to study the 
effect of glycol on the relative volatility of ethyl alcohol. The results of the investigations indicate 


the disappearance of the binary azeotrope between ethyl alcohol and water at the concentrations 


of the solvent studied in this work. 


Résumé-—Les données d’équilibre liquide—vapeur pour le systéme éthanol—eau sont déterminées 


& une pression de une atmosphére en présence d’éthyléne glycol en utilisant l'appareil d’Othmer 
de fagon a étudier l'effet du glycol sur la volatilité relative de lalcool éthylique. Les résutlats des 
recherches montrent la disparition de l'azéotrope binaire entre alcool éthylique et eau aux 


Zusammenfassung 


concentrations de solvant étudiées dans ce travail. 


Dampf—Flissigkeits-Gleichgewichte fiir das System Aethanol Wasser 


bei 1 atm Druck und in Gegenwart von Aethylenglykol wurden in einer Othmer-Apparatur 


bestimmt, un den Einfluss des Glykols auf die relative Fliichtigkeit von Aethanol zu studieren. 


Die Ergebnisse zeigen das Verschwinden des biniiren Azeotrops zwischen Aethanol und Wasser 


im Bereich der verwendeten Konzentrationen des Losemittels an. 


INTRODUCTION 


Oxpinary distillation of ethyl alcohol yields a 
product containing nearly 10-6 mole °,, of water 
due to the formation of an azeotrope. Further 
dehydration of ethyl aleohol can be carried out 
only in the presence of an added component 
using the principle of azeotropic distillation or 
extractive distillation. Benzene similar 
hydrocarbons [1], trichloroethylene [2], ether 
[3] ete., which form binary or ternary azeotropes 
with ethanol-water system have been used as 
entrainers. In the field of extractive distillation 
high boiling liquids such as glycerol, glycols ete., 
have been recommended by Maritier [4] as 
early as 1936. Crperquist [5] reports tho 
economical working of the plants in Sweden 
using ethylene glycol as the extracting agent for 
the production of absolute alcohol. However, 
a systematic investigation of vapour-liquid 
equilibria of ethanol-water systems in the presence 
of ethylene glycol has not been reported in the 
literature. The present study was, therefore, 
undertaken to determine such data which are 
essential for the design of the extractive distilla- 


tion unit. 


Usep 


MATERIALS 


Ethyl alcohol: Dehydrated alcohol supplied by the 
Madras Medical Store was distilled over calcium turnings 
and used. (Density d3° 0-7808: Refractive index np 
13570). 

Ethylene glycol: Pure glycol supplied by E. Merck 
& Co. was distilled under a vacuum of 72 em Hg and the 
fraction collected at 120°C was used. (Density d,*° 
1/1063: Refractive index 1-4283). 

Water: Laboratory distilled water was redistilled and 


used, 


EXPERIMENTAL PROCEDURE 


Vapour liquid equilibrium data were determined 
using a modified form of Othmer’s still described 
by Orumer et al. [6] with slight modification to 
improve the accuracy of the readings as evaluated 
by a series of runs with an ethyl alcohol—water 
system. The upper portion of the boiling chamber 
and the vapour outlet were provided with 
nichrome coils embedded in a coating of asbestos 
magnesia to maintain the temperature of the 
walls slightly higher than the boiling point to 
avoid errors due to condensation. 

The charge for the experimental runs was 
prepared by mixing weighed amounts of the three 
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components, so that, while the solvent concentra- 
tion was kept a constant, the other two com- 
ponents were varied, A known amount of a 
mixture corresponding to the approximate vapour 
composition was added to the charge initially 
and the three way stopcock was adjusted to 
build up this much of condensate hold-up during 
the runs. The steady conditions were reached 
after about 2-3 hr as indicated by the constant 
volume of the hold-up and the constant boiling 
temperature of the liquid. The system was 
allowed to maintain this equilibrium state for 
about 1 hr, before samples were taken in stoppered 
iodine flasks which were immediately cooled. 
The analysis of the ternary mixtures was 
based on the measurements of physical properties 
such as density and viscosity. Synthetic mixtures 
of the three components were prepared by 
accurate weighing and the physical properties 
were measured at 30°C in a thermostat. The 
readings were plotted to yield two families of 
curves for density and viscosity respectively and 
these graphs were used for analysing the equili- 
brium samples by interpolation. The small 
amounts of ethylene glycol present in the vapour 


samples were analysed by titration following the 
method described by Curme and Jounsron [7]. 


Discussion OF RESULTS 

Table 1 gives the experimental vapour liquid 
equilibrium data of ethanol-water on a glycol 
free basis at nearly constant concentrations of 
glycol in the liquid, along with the boiling points 
of the ternary mixtures. Fig. 1 shows the boiling 
point diagram for the ethanol-water system 
containing 50 per cent, 75 per cent and 90 per cent 
by weight of ethy lene glycol in the overall mixture. 
It can be observed that the boiling point curve 
tends to become flat at the high alcohol region 
for 50 per cent glycol concentration, while 
increased solvent concentration makes the curve 
steeper. The differences in boiling points between 
the water end and the alcohol end of the boiling 
point curves are 17-6, 16-2 and 15-5 °C at 50 per 
cent, 75 per cent and 90 per cent of glycol by 
weight respectively, whereas it is 21-6 °C for the 
binary ethanol-water system. Fenske et al. 
[8] have reported that the boiling point difference 


Fic. 1. Boiling point —- composition diagram. 
0 


Glycol by wt. 4 7% 
0% \ 90% 


increased with higher solvent concentration for 
the toluene-methyleyclohexane system with aniline 
as the solvent. 

Fig. 2 shows the equilibrium data for the 
ethanol-water system plotted on glycol free basis 
with the concentration of glycol as a parameter. 


Fic. 2. Vapour liquid equilibrium diagram. 
Glycol by wt. 75% 
A 9% 
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Vapour liquid equilibria for ethanol-water-ethylene glycol system 


Table 1. Vapour liquid equilibrium data of 


ethanol—water—ethylene — glycol at’ atmosphere 


pressure 


Temp. 
(©) 


107- 2 Ooo 
105-8 ‘7 | O 
| 
0-048 
0-060 
0-082 
107 
164 


95 
O-6 241 O-SAT 381 — 
265 
o24) $24 Fic. 3 
w 


Relative volatility — composition diagram. 


Glycol by wt. 0° 75 


50% 90% 


on a solvent free basis as defined by the following 


5°: 
114-0 ow equation 
24 5 i iw 
11 rl u ¥,/(1 — 
110-7 0-068 270 5-06 - 
108-5 3-6 3-5 O-104 O-361 


104-4 0 2-7 27 0-585 3-79 were determined. Fig. 3 shows the relative 

103-4 is 2-4 0-392 0-675 3-22 volatility of aleohol with respect to water plotted 
12-4 6 22 05M) 784 310 against the concentration of ethyl alcohol in the 
102-2 22 0-067 O-S45 2-72 

liquid for different solvent concentrations as 
1OL-8 3-7 24 0-939 0-991 v.15 parameter. It may be observed that the relative 
136-3 0-0 volatility is enhanced by glycol above a critical 
17-6 0-027 0-102 concentration of ethanol viz. about 24 mole °, 
Lit 88-2 14s O-047 0-166 $06 on a glycol free basis, which is common for all 
31-3 7-8 20 0-273 3-83 
‘ the three concentrations of glycol studied and is 


depressed below this critical value of alcohol. 


NOTATION 


relative volatility of alcohol to water based on 


mole fraction on glycol free basis 


glycol 


mole fraction of alcohol in the liquid on a 


free basis 


It can be observed that the effect of solvent is to X, = weight per cent of glycol in the liquid 

make the system behave more ideally. The 4) mole fraction of alcohol in the vapour on a glycol 
azeotrope of alcohol and water is broken even by free basis 

the presence of 50 per cent glycol in the solution. Y, = weight per cent of glycol in the vapour 


The relative volatility values of alcohol based T = boiling point “( 


| | | 
| 
0 
4 6} + 4 
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Notes on the diffusion-type model for longitudinal mixing in flow 


(Concerning papers by O. Levensrre., W. K. Smrru, E. Van per Laan and R. Arts) 


Anis has recently written a letter [1) increasing the value 1 dRiz.-) 1 dR, (z.*) 
of Va~ per LAAan’s extension [2] of the methods used by 
Pe dz Pe, dz 
Levenspret, and [3] to determine dispersion 
coetlicients for longitudinal mixing. In it, he showed that R (2) Ry concentration measured at = z (3) 
the impossible experimental task of injecting a perfect Van per Laan’s nomenclature and equation numbers 


delta function of tracer into the flowing system is un- are used except for (3'). 


necessary if the tracer concentration is measured at two The last boundarv condition is due to the fact that 


points in the system. As Ants mentions, the differential 9 as well as at 2, 


and the injection is at 2 ty. Upon checking the solution 


of these equations as given in Arts’ letter, it was found 


measurements are now being taken at 2 
equations to he solved are 


VOL. (2 . 1 > R a a ; : (1b) that two errors existed: (a) a misprint of the first terms 
L2 Pe ; in the numerator and denominator which should have 
1960 exponents (y Pe and (q 4) 1) re- 
> ) 1 >» spectively, and (b) a more serious error in the constant 

> Pe, “1 (Ie) coeflicients of the exponential terms ((1 q — q) and, 

' (1 —q — q)). The given solution is not consistent with 

Ry | 0 ) finite (3b) the boundary conditions equations (3f) and (3d), as may 


be readily shown by differentiating R and R, and sub- 


R (2,7) Ry, (2,7) (3f) stituting into (3f) and (3d). 


The correct solutions are obtained by solving equation (1b), (le), (3b), (3f), (3d) and (3°) by the Laplace trans- 


formation, giving 


(¢ — %) exp (4 Pe (2, 2) + exp }) Pe (2, } 
\(q — exp }) Pe (2, 29) + exp }) Pe (2, 
R 2 R, exp (YW }) Pe, =,) — (b) 
(y EXP | (y Pe (2, 9) (4 exp (4 Pe (2, =9) 


These equations may be compared with those of Arts in order to see the differences. At the second measuring 


point 2 (2 = =). the solution is found by substituting = Zp, into R (equation a). 
R R | (q — exp Pe (4 exp Pe | (c) 
where 
The following procedure will be used to find the first R, = 1 — we P + 4 eg? + wg?) p? + -- (d) 
and second moments of 2. As noted by Ants, if the 
above expression is expanded in a Maclaurin series, the Hom P (o,, Hm?) (e) 
coefficients of p™ in the expansion are 
(-1)" The expression in the braces in the formula for R, may 
(nth moment about the origin). be written as: 
If is defined as the second moment about the mean, 1 up (f) 
then the second moment about the origin is (o* p*) [4]. 
Therefore, the following expansions can be made : m 1 Pe" (1 8) 1 exp (— Pe)] exp — Pe ¢ (g) 
2 1-8 
o 3pexp (— 2 Pe [exp(— 2 Pe) 1 8) exp (— Pe [exp (— Pe) 1} 
Pe? 


t+ Peexp(— Pe (— Pe) tj} (h) 


& 
i 
= 
5 
= 


Letters to the Editors 


Therefore the values of » and o* are different than those Since the various dimensionless variables in VAN DER 


of Ants due to the error in the constants as mentioned LAAN’sS nomenclature involve the distance between 


above. Then, as shown by Arts: injection and measurement, comparing the solutions for 


R P 2 jig?) two different measurement points involves reverting to 


up (o? y2) p? absolute (dimensional) variables. 


1 — (pg + + Of + + For variables based on a length 


4 4 


where p and o? are given by (g) and (h) the absolute mean and variance are related to the dimen- 


siontless ones by 
These values of » and o- should check with those of 


Va~ per Laan for the special case of a perfect delta 


2 2 
function input. To show that this is true for my solution, Ka My % =) or (1) 
“a “ 
the following nomenclature will be used : 
ata rj ro. & perfect delta function is injected’ Treating the means first : 
ata the concentration is measured and 
Ma Moa Me 
at a _ the concentration is measured and has “ “ 
mean and varhanece 
where a4 distance from entrance of finite tube . (“= "), (m) 
ma 
(VaN per LAAN). 
He ~ B) {exp D) (4, — 4,,) exp (a, (nm) 


From Va~ per Laan’s Table T: 


fexp D) (a —exp Dy ta 
i Pi 1 I 1 


m 


xpexp | (nu D) a; (1 Byexp (u D) (a, (p) 


Therefore, jag Sheng Hoa. &S predicted by equation (j). Similarly, the variance, equation (k), may be proven to 


be correct. 


Further extensions of these methods are forthcoming. 


K. B. Biscnor 
Illinois Institute of Technology, 
Chicago, Illinois 
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The Structure and Properties of Porous Materials. 
Edited by D. H. Evererr and F.S. Stone. Butterworths, 
London 1958 and Academic Press, New York 1958. 
xis 380 pp. $12.00, 


Tue present state of knowledge on surface area of porous 
masses forms one of the main groups of subjects discussed 
at the 10th Meeting of the Colston Research Society of 
Bristol, England, and the book under review contains 
the proceedings of this symposium. This was a gathering 
of 80) surface chemists and thermodynamicists from 
eleven countries who deliberated on eighteen papers, the 
discussion takiag up almost half the book and containing 
much original research work apart from the contents of 
the papers. The structure and particularly the surface 
characteristics of crystals, gels and carbons, and the 
relevant properties of models of porous systems figure 
very largely as substance or background to most of the 
papers even when the swelling properties of wool and 
nvlon or the weathering properties of building stone are 
discussed. Only few papers contain direct application 
to reaction kinetics of flow systems or to drying ; one 
entitled Porous structure of paper deals with filtration, 
and another entitled The influence of surface charges 
deals with ion exchange. The major contributions both 
in substance of the papers and in the discussion came from 
Professor R. M. 
Professor J. H. de Borr, Technical University, Delft ; 
Professor J. R. Dacry, Royal Military College, Kingston, 
Ontario ; Professor D. H. Evenrerr, University of Bristol ; 
Dr. B. Imecik, Sorbonne, Paris ; Professor A. V. Kiseiev, 
K. Wynne- 


Jones, King’s College, Newcastle upon Tyne. 


BARRER, Imperial College, London ; 


University of Moscow and Professor W. F. 


Relief from the bewildering abstractions of hopping, slid- 
ing and sorbed molecules, of ink bottle pores and molecular 
sieves is provided by a most fitting and inspiring final 
paper by Professor Preston, Leeds University, on some 
properties of living matter, viz. the movement of sap and 
carbohydrates in large trees. 

P. ErseNKLAM 


J. S. Rowiwson: Liquids and Liquid Mixtures. 
Butterworths, London 1959. 360 pp. 75s. 


Tur title of this book is somewhat misleading as it deals 
only with equilibrium properties. The systems considered 
are generally simple condensed gases and organic liquids. 
Metals, electrolytes and polymers are not discussed. The 
subject matter can be roughly divided into three main 
statistical mechanics and 


themes, thermodynamics, 


experimental methods and results. 


The first main section deals with thermodynamic 
properties and the critical state of pure liquids. Together 
with relevant thermodynamics, vapour pressure, thermal 
expansion, isothermal compressibility, thermal pressure, 
adiabatic properties, heat capacity and residual and 
configurational properties are discussed. Throughout, 
close contact with experiment is maintained, not in the 
form of detailed descriptions of apparatus, but rather in 
the mention of important points for the experimentalist. 
Useful tables of the thermodynamic properties of a number 
of the simpler liquids are given. 

In the following two chapters on the properties of simple 
and complex liquids at low pressures the basic thermodyna- 
mies of ideal, non-ideal and partially miscible liquids is 
treated. The experimental data for approximately one 
hundred binary systems are discussed ranging from 
mixtures of simple condensed gases to aqueous mixtures, 
The thermodynamic treatment of consolute temperatures 
which is based on a Taylor expansion of the free energy 
about its value at the critical point is similar to that 
used for discussing critical phenomena in pure liquids. 
Testing of thermodynamic data for consistency is stressed. 

The properties of liquid mixtures at high pressures is 
treated under azeotropy, latent heats of binary mixtures 
and critical phenomena in mixtures of simple and complex 
liquids. Much of the older work in this field is usefully 
summarized here. 

The last main section dealing with intermolecular forces 
and the statistical mechanics of fluids and their mixtures 
ean be read independently of the previous chapters. It 
provides an up-to-date account of the various attempts 
which are being made to calculate the equilibrium pro- 
perties of liquids. A not unwelcome note of realism can 
be detected by the reviewer in this section when the 
author, facing realities, holds out that the modified cell 
theories provide more immediate hope for progress than 
the distribution function approach which is still anchored 
to the superposition approximation. The book is well 
produced, liberally illustrated and has an_ extensive 
bibliography. However the use of lower case letters, 
inconsistently applied to the titles of journals listed in the 
references, is in opposition to the currently accepted format 
and has nothing to recommend it. The number of printing 
errors is a minimum. 

As a source of information on the experimental equili- 
brium properties of liquids and liquid mixtures and an 
up-to-date account of the associated thermodynamic 
and molecular aspects of the subject, this book can be 
confidently recommended to engineers. The price is 
unfortunately rather on the high side. 
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Isotopentechnik in Amerika. Herausgegeben vom 
Rationalisierungs-Kuratorium der Deutschen Wirtschaft. 
Auslandsdienst, Helft 84. Carl Hanser Verlag, Miinchen 
1959. 1381S. 17 DM. 


Dem Bericht liegt cine Reise von Vertretern der Wissen- 
schaft und Wirtschaft, der BehOrden und Gewerkschaften 
Berlins zu Grunde, in deren Verlauf Ende 1956 in den USA 
mehr als 20 Unternehmen und Institute besucht wurden 
die sich mit Isotopen- und Kernenergietechnik befassen 
An allgemeinen Problemen werden die Herstellung 
radioaktiver Isotope, Fragen der Laboreinrichtung, 
Sicherheitsvorschriften und Massnahmen, der Transport 
und die Abfallbeseitigung besprochen. \nschliessend 
werden — unterstiitzt von Abbildungen und Literatur- 
hinweisen — Anwendungen von Isotopen in den verschie- 
denen Zweigen der Industrie und Strahlungsmessgeriite 
beschrieben und in aller Kiirze einige EKinblicke in die 
Reaktortechnik gegeben Kin besonderer Abschnitt ist 
der wirtschaftlichen Bedeutung der Kerntechnik gewidmet 
Ohne sich in technische oder physikalische Einzelheiten 
zu verlieren, gibt das leicht lesbare Heft cinen Kindruck 
von dem Stand der Entwicklung und ermoglicht einen 
Uberlick fiber die Bedeutung, die Probleme und die 
Anwendungsmoglichkeiten der Isoitopen- und Kernenergie- 
technik. 
G. Hevi 


Grundlagen der Chemischen Technik (Herausgeber, 
Hl. und Fuens) Band Fucns: 
Physikalische Chemie als Einfiihrung in die 
chemische Technik. Verlag Sauverlinder, \arau und 
Frankfurt (Main) 1959. 4968S. 37 DM. 


Der zweite Band der Buchreihe Grundlagen der Chemischen 
Technik mit einem der Herausgeber als Autor behandelt 
diejenigen Kapitel der physikalischen Chemie, welche 
hiufig als Grundlagen in der chemischen Technik ge- 
braucht werden. Er tut dies in der richtigen Erkenntnis, 
dass insbesondere der Verfahrenstechniker in den physika- 
lischen Grundvorgiingen heimisch sein muss. Das Buch 
erhebt nicht den Anspruch, ein Lehrbuck oder Nachschlage- 
werk fiir physikalische Chemie zu sein, zeigt aber im 
Gegensatz zu diesen klassischen Werken in ciner Reihe von 
praktischen Beispielen, in welcher Weise die Kenntnis der 
physikalisch-chemischen Grundbegriffe unentbehrlich ist 
und wie die diesbeziiglichen Daten praktisch mithelfen, 
technische Werte zu bestimmen und damit tiefer in die 
Vorgiinge, wie sie sich in industriellen Prozessen volziehen, 
einzudringen. 

Bei der Fille des zu behandelnden Stoffes war es ohne 
Zweifel schwierig, die Grenzen fiir das Wesentliche zu 
ziehen, und einige Unterkapitel wie z.B. Explosionsgrenzen 
oder gewisse technische Abschnitte im Kapitel Phasen- 
grenzflichen mussten reichlich kurz gefasst werden. 

Das Buch enthilt erstaunlich viel Material sowohl fiir 
Studierende in Verfahrenstechnik und chemischer Techno- 


Book Reviews 


logie wie auch fiir Chemiker und Ingenieure, welche in der 
Praxis stehen und in nicht zu langer Zeit ihr theoretisches 
Wissen wieder auffrischen und ergiinzen wollen. Das 
Werk enthilt viele persOnliche Erfahrungen des Autors ; 
die Art der Darstellung und der gewihlte Stil machen das 
Studium des Buches recht leicht. 

Dem Verlag Sauerlinder gebihrt Anerkennung fiir 


die gute Ausfiihrung des Buches bei missigem Preis. 


\. GUYER 


Crype Orr Jr. and J. M. Fine Particle 
Measurement. Macmillan, New York, 1959. ‘xiv 353 
Pp. £3 libs iid. 


Tue late Professor DALLAVALLE has been a pioneer in a 
field which has now come to be known as Powder 
Technology and his well-known book Micromeritics 
remained a guide on the behaviour and characteristics of 
small particles since its first edition appeared sixteen 
vears ago. He and his co-worker at the Georgia Institute 
of Technology in Atlanta have now enlarged the section 
on particle measurement and deal with Size, Surface and 
Pore Volume. 

The book consists of a collection of numerous methods 
which have been used or suggested, ranging from micro- 
scopy and sieving to techniques making use of radiation 
scattering and transmission. No definition of “ fine 
particles” is given apart from the statement that half 
the universe consists of them, and the electron microscope 
(lower limit 0-001 «) and adsorption methods are treated 
together with sieving methods so that the size ratio under 
consideration is of the order of 1: 107. No wonder that 
within 300 odd pages only a catalogue of methods 
could be given with little attention to the relative im- 
portance or merits of each. Thus we find that elutriation 
is dismissed in 500 words, capillary pressure methods of 
surface determination in half this amount, statistical 
methods are not mentioned at all, and some methods are 
treated only in the from of cursory summaries. 

The great virtue of the presentation lies in the bringing 
together of a mass of information and treating it up to a 
level which is fully adequate for a general understanding 
of the principles and an appreciation of the use of each 
method. Some 700 references are provided as sources of 
detailed information. Although, in general, these are 
adequate, one misses some people well known on this side 
of the Atlantic, foremost being Dr. H. Hrywoopn, an 
authority on size measurement of small particles since the 
early thirties. There is evidently great uncertainty 
regarding the accuracy of interpretations of measurements 
particularly of surface area because of the lack of precise 
definitions, and hence the collection of typical results 
and comparison measurements in the appendix to the 
book is valuable. 


EisenKLAM 
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Trockungstechnik O. und 
K. Band I.—Trockner und Trockuungs- 
verfahren (K. Springer-Verlag, Berlin 1959. 
S88p. 69.00 DM. 


Der zweite Band der Trocknungstechnik behandelt die 
technischen Verfahren und Apparate zur Trocknung 
bestimmter Giiter unter Beriicksichtigung der im ersten 
Band beschriebenen wissenschaftlichen Grundlagen des 
Trocknungsvorganges. 

Der umfangreiche Stoff, nach dem  Gesichtspunkt 


gleicher physikalischer Gesetzmissigkeit gegliedert, wird 


dem Leser in tibersichtlicher Weise geboten. Dabei wird 
jede Trocknerbauart sowohl aus der mechanischen, 
thermischen und aerodynamischen Wirkungsweise des 
Trockners als auch aus den technologischen Kigenschaften 
des darin behandelten Gutes heraus verstindlich gemacht. 

Das Buch zeigt am Anfang, wie man Trockner planen 
und entwerfen kann, und gibt Hinweise fiir die Umrechnuny 
von Versuchsergebnissen auf andere Temperaturen, rela- 
tive Feuchten, Schichtdicken usw. Ein weiteres Kapitel 
ist dem Wiarmebedarf des Trockners und der Frage der 
Wirmeriickgewinnung gewidmet. Im Abschnitt “Férdern 
und Fithren der Luft © wird durch zahlreiche, anschauliche 
StrOmungsaufnahmen am Wasserkanal-Modell die stré- 


mungsygerechte Gestaltung von Trocknern gezeigt. Ein 


umfangreiches Kapitel tiber Konvektionstrockner behan- 
delt die verschiedenen Mdglichkeiten der Um — bzw. 
DurchstrOmung des Gutes und beschreibt neben einer 
Vielzahl spezieller Bauformen auch die Trocknung be- 
stimmter Giter, so z.B. die Holz-, Leder-, Lack-, Textil-. 
Getreide-, Lebensmittel- und Futtermitteltrocknung. Die 
Verfahren der Strahlungstrocknung und der Trocknung im 
elektrischen Wechselfeld werden ausfiihrlich nach physika- 
lischer Wirkungsweise und technischer Ausfiihrungsform 
diskutiert. Das Buch befasst sich ferner mit der Technik 
der Kontakttrocknung und der Vakuumtrocknung ein- 
schliesslich der Gefriertrocknung. Ein besonderes Kapitel 
zeigt Moglichkeiten zur Luftentfeuchtung und zur Wieder- 
gewinnung von Losungsmitteln, die beim Trocknungs- 
vorgang aus dem Gut verdunsten. 

Das Werk gibt, zusammen mit dem ersten Band. 
einen umfassenden Cberblick alle mit der Trocknunyg 
zusammenhingenden Fragen und bietet sowohl dem 
Hersteller als auch dem Anwender von Trocknungsanlagen 
wertvolle Arbeitsunterlagen. 

Das ausfiihrliche Inhalts- und Sachverzeichnis und 
die Schrifttumsangaben am Schluss jeden Kapitels tragen 
zur Ubersichtlichkeit des Buches wesentlich bei, das 
wegen der leicht verstindlichen Fassung auch dem 
Studierenden sehr empfohlen werden kann. 

H. Kunz 
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SELECTION OF CURRENT PAPERS OF INTEREST TO 
CHEMICAL ENGINEERS 


W. M. Jones and P. J. Isaac: The flow of carbon dioxide and nitrogen at high pressures through porous plugs of lamp- 
black. Trans. Faraday Soc, 1959 55 19047-1958. 


G. 1. Tayvtor: The dynamics of thin sheets of fluids—-I. Water bells (shape produced by obstruction in a jet), IL, Waves 
on fluid sheets. IIL. Disintegration of fluid sheets. Proc. Roy. Soc. 1959 A253 289-321. 


G. A. and R. C, LI. Jenkins Sedimentation and electrophoresis of interacting substances IL. Asymptotic 
boundary shape for two substances reacting irreversibly. Proc, Roy. Soc, 1959 \253 420-437 (See also G, A. GILBERT : 
Proc. Roy. Soc. 1959 A250 377-388. 


Lu Ting: On the mixing of two parallel streams. J. Math. Phys. 1950 38 153-165 
s. W. Yuan: Turbulent flow in channels with porous walls. J. Math. Phys. 1959 38 166-171 
G. D. Gorvon : Mechanism of speed and break-up of drops (in an air stream). J. Appl. Phys, 1959 30 1759-1761, 


I. S. Bsorkiunp and W. M. Kays: Heat transfer between concentric rotating cvlinders. J. Meat Transfer (Trans 
imer. Soc. Mech. Engrs. Series C) 1950 81 175-186. 


J. KE. Lay: An experimental and analytical study of vortex flow temperature separation by superposition of spiral 


and axial flows (Parts 1 and 2). J. /leat Transfer (Trans. Amer. Soc. Mech. Engrs. Series C) 1959 81 202-222. 


kK. M. Sparrow and J. L. Greoa: Heat transfer from a rotating disk to Nuids of anv Prandtl number. J. Meat Transfer 
(Trans. Amer. Soc. Mech. Engrs. Series ©) 1959 81 249-251 


V.S. Anpact and J. A. CLark : Dynamic response of heat exchangers having internal heat sources—II1. J. Meat Transfer 
(Trans. Amer. Soc. Mech. Engrs. Series ©) 1959 81 253-266 (See also Trans. Amer. Soc. Mech. Engrs. 1958 80 612-024, 
Ibid. 1958 80 625-634). 


Hl. Wour: Heating and cooling air and carbon dioxide in the thermal entrance region of a circular duct with large gas 
to wall temperature differences. J. lleat Transfer (Trans. Amer. Soe. Vech. Engrs. Series C) 1959 81 267-279. 


KR. Steere. and KE. M. Sparrow: Turbulent flow in a circular tube with arbitrary internal heat sources and wall heat 
transfer. J. Heat Transfer (Trans Imer. Soc. Mech. Engrs. Series CC) 1959 81 Zso-200. 


kK. M. Sparrow and J. L. Greec : Laminar condensation heat transfer on a horizontal evlinder, J. Meat Transfer (Trans. 


Amer. Soc. Mech. Engrs. Series ©) 1959 81 291-296. 


kK. L. Lustenaper, R. Ricwrer and F. J. Neccesaver: The use of thin films for increasing evaporation and condensa- 


tion rates in process equipment. J. Heat Transfer (Trans. Amer. Soc, Mech. Engrs. Series ©) 1959 81 297-307, 


KR. L. Cuampers and FE. V. Somers: Radiation fin efficiency for one-dimensional heat flow in a circular fin. J. Heat 


Transfer (Trans. Amer, Soc. Mech. Engrs. Series C) 1959 81 327-329. 
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SELECTION OF CURRENT SOVIET PAPERS OF INTEREST TO 
CHEMICAL ENGINEERS* 


EK. K. Sirpe and P. G. Romankov: Investigation into the process of steam distillation. Zh. prikl. Khim. 1959 32 


2197-2207. 


Investigation of the rate of absorption of 


S. N. Ganz, S. B. Lemovicn, N. A. and M. A. : 
carbon dioxide by monoethanolamine in a horizontal mechanical absorber. Zh. prikl. Khim. 1959 32 2207-2210. 


A. A. Noskov, G. V. Burova and P. Feipresn: Plate efficiency of a single sieve plate in rectification. Zh. prikl. Khim 


1959 32 2211-2218. 


V. B. Sarkrrs, D. T. Traper and I. P. MukHLenov: Heat transfer from a suspended catalyst bed to a heat transfer 


surface. Zh. prikl. Khim. 1959 32 2218-2225. 
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1960 


B. M. BoGos avs, B. M. Yavorski and A. D. Virntk: On application of thermal diffusion to purification of dyestuffs 


Z. prikl. Khim, 1959 32 2225-2229. 


M. V. Pourakoyv, V. V. Suarya and Z. Z. Visorski: Investigation of catalytic conversion of methanol into formal- 


dehyde in a fluidized catalyst bed. Zh. prikl. Khim. 1959 32 2275-2283. 


D. P. Dosicuin: Diffusion of water vapour and flow of air in porous glass. Zh. prikl. Khim. 1959 32 2336-2339. 


Khim. Nauka i Prom. 1959 4 554-565. 


P. A. Resinper: Surface-active substances and their application. 


Physico-chemical principles of wetting and of detergent action of surface-active substances. Ahim 


\. B. TaupMAN : 
Nauka i Prom, 1959 4 566-573. 


Methods of synthesis of surface-active substances. Khim. Nauka i Prom, 1959 4 


A. D. Perrov and G. I. : 


573-585. 


F. V. Nevouin: Synthetic detergents and their compositions. Khim. Nauka i Prom. 1959 4 586-592. 


R. A. Patina: Surface-active substances on ethylene-oxide basis. Ahim. Nauka i Prom. 1959 4 592-598. 


A. K. Livsarrs: Flotation reagents in current use, their production and application. Khim. Nauka i Prom. 1959 4 


622 G28. 


M. M. Ecorov, V. F. Kisecev and K. G. Krasttnixov: Adsorption capacity of unit surface of quartz. Zh. fiz. Khim. 
1950 33 2141-2144. 


M. G. KaGaner: New method of determination of specific surface area of adsorbents and other finely-dispersed sub- 


stances. Zh. fiz. Khim. 1959 33 2202-2210. 


Kinetics of attachment of bubbles to mineral particles during 


B. V. Dervacury, 8S. 8S. and V. A. LisicHenko : 
flotation. Zh. fiz. Khim. 1959 33 2280-2287. 


Equations of extraction 


Theory of diffusional extraction of substances from porous bodies—TI. 


G. A, AKSELROD : 
kinetics. Zh. fiz. Khim. 1959 33 2316-2324. 


M. M. Dusinin and E. G. Zaukovskaya: On adsorption properties of carbon adsorbents. Study of adsorption pro- 
Im. Akad, Nauk SSSR, Otd. khim. Nauk 1959 


perties of active charcoals with developed transitional porosity. 


(10) 1705-1715. 


Estimation of influence of free convection on turbulent flow. Zh. tekh. Fiz. 1959 29 1390-1392. 


E. YANTOVSKI: 
Dokl Akad. Nauk SSSR 1959 


G. D. RosenserG: Experimental investigation of unsteady flow of viscoplastic fluids. 
129 56-58. 


*To assist readers, translations of any article appearing in the above list can be obtained at a reasonable charge. All orders should be addressed 
to the Administrative Secretary of the Pergamon Institute at either Headington Hill Hall, Oxford or 122 East 55th Street, New York 22 which- 


ever is more convenient. 
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Selection of Current S viet Papers of Interest to Chemical Engineers 
F. Kazantsev: Motion of air bubbles in water under the action of Bjerknes forces in acoustic fields. Dokl. Akad. 
Nauk SSSR 1959 129 64-467. 


’. Yu Onvov and N. M. ZHavornonkovy Effect of ultrasound on the absorption of carbon dioxide by water. Dokl. 
thad. Nauk SSSR 1959 129 161-164 


P. Timorreyv and |. T. ErxasnKko: Investigation of the structure of sorbents by kinetic methods. Dokl, Akad. Nauk 
SSSR 1959 129 384-386 


L.. KuTMAN : Effect of speed of rotation and of shape of armour on performance of a drum-type ball mill. Teploener- 


getika 1950 © (11) 24-28. 


F. Drryakts and L. N. Barrxeva: Generalizacion of measurements of drop sizes with the aid of dimensionless criteria 
in the case of atomization of liquids by a centrifugal atomizer. Teploenergetika 1959 6 (11) 33-36. 


V. Sono ov and V. P. Osoxrs \brasivity of coal and the resistance of metals to wear in milling. Teploenergetika 
1959 © (11) 37-41 


P. Vokovovicn and V. V. Avrunin: Experimental investigation into the P-V-T relationship of carbon dioxide 


Teploenergetika 1959 6 (11) 58 65 


ERRATA 


SELECTION OF CURRENT Sovier Parers or Inrerest ro Cuemicat Engineers. Chem. Engng. Sci. 1959 11 222 


The last reference on p. 223 should read: S. S. Kuraretapze and V. N. Mosxvicneva: On the relation between 
the hydrodynamics of a two-component layer and the theory of crisis in the mechanism of boiling. Zh. tekh. Fiz. 1959 
29 1133-1139 


The first two references on p. 224 should read: P. N. Kupansxi: On the coagulation of aerosols by free streams 
of gas. Zh. tekh. Fiz, 1959 29 1140-1141. V. V. Racwinskaya: On the problem of hydrodynamic properties of ion 
exchange resins. Zh. tehh. Fiz. 1959 29 1159-1161. 
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